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Abstract 


A  technique  for  recording- and  visualizing  acoustic  wavefronts  over  a  finite  field 
of  view  and  with  a  resolution  on  the  order  of  Angstroms  has  been  developed. 
Multiple  exposure  holograms  of  an  aluminum  test  black  were  used  to  record  the 
surface  displacements  resulting  from  the  detonation  of  a  small  explosive  charge 
placed  on  the  object  surface.  The  primary  wavefronts  resulting  from  the  explosive 
charge  were  sufficiently  large  as  to  be  visible  upon  reconstruction  of  the  holograph¬ 
ic  interferograrns.  Tn  addition,  details  of  the  surface  wavefronts  which  could  not 
be  deduced  using  conventional  holographic  interferometry  were  measured  and 
displayed  by  modifying  the  pulsed  recording  geometry  and  subsequently  applying 
heterodyne  interferometric  analysis  to  the  holograms  to  interpolate  between  the 
observed  interferometric  fringes.  Although  theory  predicts  that  still  greater  resolu¬ 
tion  is  possible,  fringe  interpolation  to  1/900  of  a  fringe  (about  3a  displacement) 
was  demonstrated.  Rased  on  these  experimental  results  and  the  mathematical  rela¬ 
tionships  which  have  been  developed  to  describe  system  performance,  it  appears 
likely  that  a  practical  system  for  studying  sub- Angstrom  transient  phenomena  may 
be  realized  using  heterodyne  holographic  interferometry. 
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INTRODUCTION 


To  dale,  conventional  techniques  for  measuring  and  studying  transient  surface 
wave  phenonena  have  been  limited  to  transducer  techniques  which  provide  a  time 
record  of  surface  displacement  at  a  single  position  on  an  object  surface.  Piezoelec¬ 
tric,  capacitive,  magnetic,  and  single  point  interferometric  techniques  all  fall  into 
this  category.  For  many  testing  and  measurement  applications,  these  techniques 
are  quite  adequate.  However,  some  studies  are  difficult  if  not  impossible  to  per¬ 
form  with  such  transducer  systems.  If  one  wishes  to  investigate  the  interaction  of 
surface  waves  with  a  surface  defect,  for  example,  a  single  point  transducer  must  be 
moved  around  the  defect  area  and  the  experiment  repeated  for  each  new  transduc¬ 
er  location  in  order  to  map  out  changes  in  the  amplitude  and  shape  of  the  surface 
wavefront  as  it  is  diffracted  and  attenuated  by  the  defect  For  some  transducer 
types  the  effects  of  mounting  the  transducer  on  the  surface  must  also  be  taken  into 
account.  Further,  the  acoustic  wave  source  must  be  very  reproducible  so  that  one 
is  assured  that  for  each  new  placement  of  the  transducer,  the  wave  which  is  pro¬ 
pagated  is  identical  to  the  one  before  it  Alternatively,  one  can  imagine  an  experi- 
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roent  where  an  array  of  noncontact  ing  transducers  are  positioned  on  an  object  sur¬ 
face  so  that  surface  displacement  over  a  fixed  field  may  be  recorded  as  a  function 
of  time  and  position  for  a  single  acoustic  event.  The  cost  of  implementing  such  a 
system  and  restrictions  on  object  surface  geometry  would  seem  to  preclude  the  use 
of  such  a  system  as  a  practical  measurement  technique. 

Holographic  techniques,  on  the  other  hand,  provide  a  means  to  interfemmetri- 
cally  compare  changes  in  surface  topography  of  objects  of  arbitrary  shape  over  a 
finite  field  of  view.  Time-averaged  holographic  interferometry  has  been  used  to 
study  the  rather  large  surface  displacements  associated  with  mechanical  resonances. 
Transient  surface  displacements  associated:  with  impact  testing  and  high  energy  ul¬ 
trasonic  testing  have  been  recorded  using  double  exposure  pulsed  holography.  Tn 
both  of  these  holographic  techniques,  interferometric  fringes  are  observed  in  the 
reconstructed  holographic  image  wherever  the  surface  is  displaced  by  an  amount 
greater  than  one  half  the  wavelength  of  the  illuminating  laser  light.  Thus  displace¬ 
ments  less  than  about  250nm  (2500 A)  do  not  give  rise  to  holographic  fringes. 
Consequently  sub-fringe  wavefronts  or  sub-fringe  detail  of  large  wavefronts  cannot 
be  discerned  by  either  of  these  holographic  techniques.  However,  by  modifying 
the  holographic  recording  geometry  and  applying  heterodyne  analysis  techniques  to 
the  recorded  hologram,  full  field  measurement  and  visualization  of  sub-fringe  sur¬ 
face  displacements  such  as  those  associated  with  traveling  surface  acoustic  wave 
may  be  achieved. 

Full  field  recording  is  achieved  using  pulsed  holography  to  'freeze"  surface  mo¬ 
tion  at  a  particular  instant  in  time.  Using  double  exposure  holography,  the  surface 


topography  of  a  test  specimen  at  different  instants  in  time  may  be  interferometri- 
c&ily  compared.  The  double  exposure  hologram  records  two  virtual  images  of  the 
test  specimen  which,  when  reconstructed  by  reillumination  of  the  hologram  with 
its  reference  beam,  appear  in  space  at  the  original  object  location  relative  to  the 
hologram  plate.  These  two  virtual  images  may  be  used  as  virtual  objects  for  an  in¬ 
terferometer  if  reimaged  on  an  observation  plane.  The  interference  pattern  pro¬ 
duced  at  the  observation  plane  is  identical  to  that  which  would  be  observed  if  the 
light  reflecting  from  the  surfaces  of  two  nearly  identical  real  objects  were  combined 
by  substituting  the  object  surfaces  for  the  minors  of  a  Michelson  interferometer. 

This  Michelson  interferometer  analogy  for  holographic  interferometry  proves 
to  be  useful  in  predicting  the  effects  of  surface  displacements  on  the  observed 
holographic  interference  pattern  In  fact  the  results  obtained  using  this  analogy  are 
exact  for  the  case  where  both  the  viewing  angle  and  the  angle  of  object  illumina¬ 
tion  are  normal  to  the  object  surface  when  the  hologram  is  exposed  and  where 
plane  wave  illumination  is  used.  Assuming  then  that  this  is  indeed  the  case  one 
can  see  that  local  changes  in  surface  displacement  alter  the  optical  pathlength  from 
the  object  surface  to  the  point  of  observation  If  at  a  single  point  on  the  object 
there  is  no  displacement  between  the  first  and  second  holographic  exposures,  then 
the  light  reflecting  from  corresponding  points  of  the  two  virtual  objects  travels  ex¬ 
actly  the  same  distance  to  the  observation  plane.  Consequently,  the  electric  fields 
of  the  two  reflected  rays  will  be  in  phase  and  interfere  constructively  to  give  a 
bright  point  in  the  observation  plane.  If,  on  the  other  hand,  at  a  different  point  on 
the  object,  some  displacement  occurred  between  exposures,  then  the  light 
reflecting  from  corresponding  points  on  the  two  virtual  objects  will  not  be  in  phase 


at  the  observation  plane.  If  the  displacement  which  occurred  between  exposures 
was  sufficient  to  cause  an  optical  pathlength  change  of  one  half  the  wavelength  of 
the  illuminating  light,  then  the  two  rays  from  the  corresponding  points  on  the  vir¬ 
tual  object  will  differ  in  phase  by  180*  thus  interferring  destructively  giving  a  dark 
point  in  the  observation  plane.  As  in  the  case  of  a  Michelson  interferometer, 
object  motion  of  one-half  wavelength  yields  one  whole  wavelength  of  optical  path- 
length  change.  Therefore  the  interferometric  fringes  formed  over  the  image  of  the 
object  in  the  observation  plane  appear  to  map  out  contours  of  constant  displace¬ 
ment  where  the  contour  interval  is  one  half  the  wavelength  of  the  illuminating 
light. 

With  the  establishment  of  this  half  wavelength  contour,  or  fringe,  interval  the 
sensitivity  limit  of  visually  interpreted  holographic  interferograms  is  set.  In  the  ex¬ 
periment  reported  in  this  work,  an  argon  ion  laser  with  an  output  wavelength  of 
514.5nm  was  used.  Thus  one  full  fringe  corresponds  to  a  surface  displacement  of 
237.3nm.  One  can  readily  discern  one-half  a  fringe  and  therefore  may  resolve 
down  to  128.6nm  by  simply  counting  fringes  observed  in  the  output  plane  of  the 
holographic  interferometer.  Unfortunately  surface  displacement  amplitudes  gen¬ 
erated  by  conventional  ultrasonic  transducers  may  range  from  0.01  nm  to  tens  of 
nanometers  and  are,  therefore,  well  below  the  limit  of  resolution  for  visually  inter¬ 
preted  holographic  interferometry.  In  other  words,  if  a  flat  mirror  were  used  as  a 
test  object  and  a  hologram  was  exposed  before  and  during  excitation  of  the  mir¬ 
rored  object  by  a  conventional  ultrasonic  transducer,  no  fringes  would  be  observed 
in  the  observation  plane  of  the  holographic  interferometer. 
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Sub-fringe  displacement  information,  however,  is  recorded  by  the  hologram 
and  may  be  extracted  using  heterodyne  interferometric  techniques  to  read  out  the 
holographic  record-  Heterodyne  techniques  employ  any  of  several  methods  to  en¬ 
sure  that  light  arriving  at  the  output  plane  from  one  virtual  object,  is  at  a  slightly 
different  frequency  from  the  light  arriving  from  the  other  virtual  object.  The 
effect  of  frequency  shifting  one  of  the  light  beams  is  to  cause  the  intensity  of  the 
image  in  the  output  plane  to  vary  periodically  at  the  difference  frequency  of  the 
two  beams.  While  the  frequency  of  the  intensity  variation  is  the  same  at  every 
point  on  the  object  surface,  different  points  in  the  image  plane  may  not  fluctuate 
precisely  in  phase  with  each  other.  The  degree  of  phase  difference  in  intensity 
fluctuation  between  two  points  in  the  observation  plane  is  directly  related  to  the 
difference  in  surface  displacement  between  the  two  paints.  Thus  if  a  one- half 
fringe  difference  were  observed  between  two  points  using  conventional  holographic 
interferometry,  there  would  be  a  180*  phase  difference  in  intensity  fluctuation 
measured  between  the  same  two  points  in  the  heterodyne  case.  Since  one  fringe 
corresponds  to  a  360°  phase  difference,  then  each  degree  of  phase  difference 
corresponds  to  an  object  displacement  of  2S7.3nm  x  1/360  or  about  0.7nm  for  the 
514.5nm  line  of  the  argon  ion  laser.  Note  that  sensitivities  in  this  range  should 
permit  detection  of  displacements  from  conventional  ultrasonic  sources.  The  abili¬ 
ty  to  resolve  much  smaller  displacements  is  limited  by  factors  which  may  be  op¬ 
timized  for  a  particular  application.  These  factors  include  quantum  noise  at  the 
detectors,  speckle  noise  in  the  image,  and  environmental  effects. 
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BACKGROUND 


The  advantages  of  using  pulsed  holography  and  heterodyne  analysis  for  study¬ 
ing  traveling  surface  acoustic  waves  become  apparent  once  the  operating  principles 
are  well  understood.  To  obtain  such  an  understanding  it  is  useful  to  review  the 
principles  of  classical  interferometry  and  holography.  This  section  includes  a  dis¬ 
cussion  of  interferometry  as  a  tool  for  measuring  dynamic  surface  displacement 
The  Michel  son  interferometer  is  used  as  an  example  with  special  attention  paid  to 
noise  sources  which  limit  displacement  sensitivity  and  are  to  some  degree  common 
to  the  heterodyne  holographic  system.  These  noise  sources  include  pathlength  in¬ 
stabilities,  detector  shot  noise,  and  thermal  noise  in  amplifiers.  Heterodyne  inter¬ 
ferometry,  the  nature  of  heterodyne  signals,  and  its  effect  on  overall  signal  to 
noise  ratio  is  discussed. 

Holography  relieves  the  restrictions  imposed  by  classical  interferometry  on  the 
types  of  test  objects  that  may  be  analyzed.  Typically  one  may  inspect  flatness  of 
mirrors  or  symmetry  of  convex  or  concave  optical  surfaces  using  classical  inter¬ 
ferometers.  Using  holography,  however,  one  may  inspect  more  complex  and  even 
diffuse  surfaces.  A  hologram  may  be  used  to  reconstruct  any  arbitrary  reference 
wavefront  especially  that  reflected  by  the  test  object  at  some  earlier  point  in  time. 
This  reference  wavefront  may  then  be  combined  with  the  object  wavefront  to  pro¬ 
duce  an  interference  pattern  over  the  image  of  the  object  indicative  of  surface  dis¬ 
placements  which  may  have  occurred  subsequent  to  the  initial  holographic  expo- 


sure.  Thus,  the  entire  object  within  the  field  of  view  of  the  holographic  inter¬ 
ferometer  may  be  holographically  recorded  and  interferaraetricaily  compared.  The 
principles  of  holographic  wavefront  reconstruction,  the  use  of  pulsed  recording 
techniques  for  freezing  rapidly  occurring  displacements,  and  the  application  of 
heterodyne  techniques  to  improve  holographic  sensitivity  will  be  discussed. 


Interferometry 

The  Michelson  interferometer  shown  schematically  in  Figure  1  will  be  used  to  il¬ 
lustrate  the  interferometric  principles  applied  to  surface  displacement  measure¬ 
ments.  Tn  every  case  it  will  be  assumed  that  coherent,  single  mode,  polarized  laser 
illumination  is  used.  The  laser  beam  entering  the  interferometer  from  the  left  is 
split  by  partial  reflectance  into  two  beams  at  the  beamsplitter.  The  portion  of  the 
beam  deflected  upwards  in  the  figure  is  reflected  by  the  reference  mirror  back 
through  the  beamsplitter  and  onto  &  photodetector.  That  portion  of  the  laser  beam 
not  initially  deflected  by  the  beamsplitter  is  reflected  from  the  test  object  surface 
and  returned  to  the  splitter  where  it  is  combined  with  the  reference  beam  and 
deflected  onto  the  photodetector. 

The  object  and  reference  waves  arriving  at  the  detector  may  be  expressed  by 
their  respective  analytic  signal  representations  (Yariv,  1976)  which  is  the  harmonic 
solution  to  the  free  space  optical  wave  equation  for  the  electric  field  at  the  detector 
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where 

UR  —  reference  beam 
UQ  —  object  beam 

—  field  amplitude  of  reference  beam 
A0  —  field  amplitude  of  object  beam 
<w  —  radian  frequency  of  the  laser 
t  ”  time 

k  “  propagation  constant  2n/k 
dl  —  total  pathlength  difference 
between  the  reference 
and  object  beams. 

These  expressions  represent  the  amplitudes  of  the  two  optical  beams.  Photo- 
detectors,  however  are  square-law  devices  sensitive  to  radiant  intensity  not  ampli¬ 
tude.  Therefore  the  electrical  signal  generated  by  the  detector  will  be  proportional 
to  the  light  intensity  incident  on  the  detector. 


/DHCfo-K/*l2 
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uo  —  complex  conjugate  of  UQ 
UR*—  complex  conjugate  of  UR. 


Thus, 


—Aj+Ait+lAo  A Kco&  kdl 

-l+2-!%^cosArf/. 

A^-Ai 


If  Aq-  A^,  then 

fo-WHl-hxw  kdl).  (3) 

Since  in  both  legs  of  the  interferometer,  the  beam  is  folded  back  on  itself,  the 
pathlength  difference,^,  resulting  from  a  normal  surface  displacement,  D,  is  twice 
the  displacement  amount.  Equation  3  may  now  be  rewritten: 

Id-IAM+cos  IkD ) 

-2AK\+cos(~-D))  (4) 

It  can  be  seen  then  that  minima  in  the  light  intensity  Ming  on  the  detector  occur 

when  23—  n-^-,  n  —  1,2,3... 

4 


Although  the  Michelson  interferometer  in  this  basic  form  has  been  used  to 
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detect  acoustic  signals  ( Wellman,  1980),  several  modifications  may  be  made  to  im¬ 
prove  its  performance.  By  twice  differentiating  Equation  4  with  respect  to  D  one 
observes  that  the  greatest  change  in  intensity  for  a  given  surface  deflection  less 
than  A/2  occurs  when  D  —  (2n+l)A/8.  One  approach,  then,  to  obtaining  max¬ 
imum  sensitivity  to  surface  displacements  is  to  actively  stabilize  the  optical  path 
difference  between  the  object  and  reference  legs  using  an  electro-mechanical  feed¬ 
back  arrangement  to  drive  the  reference  mirror.  An  additional  bendit  realized  by 
using  an  electronmechanical  feedback  system  is  that  the  surface  displacement  sig¬ 
nals  of  interest  appear  as  either  control  or  error  signals  which  may  be  easily  picked 
off  of  the  feedback  control  circuit.  The  mirror  drive  control  signal  follows  and  may 
therefore  be  used  as  a  measure  of  low  frequency  surface  displacements  (Leiner, 
1978).  The  feedback  error  signal,  on  the  other  hand,  tracks  pathlength  variations 
resulting  from  acoustic  waves,  for  example,  which  occur  at  frequencies  beyond  the 
response  limit  of  the  feedback  system  (Palmer,  1977;  Palmer,  1979;  Smeets,  1977; 
Emmony,  1982).  Techniques  other  than  feedback  stabilization  have  been  used  to 
ensure  peak  sensitivity  of  the  interferometer.  In  one  system  developed  at  RCA 
labs  (Mezrich,  1974)  the  reference  mirror  is  moved  at  a  low  frequency  and  rela¬ 
tively  large  amplitude  to  assure  that  at  some  point  in  the  measurement  period  the 
interferometer  is  in  its  optimum  configuration.  This  technique  works  espesially 
well  for  measuring  continuous,  periodic  signals.  In  such  cases  the  largest  signal 
recorded  during  the  test  interval  corresponds  to  the  point  of  maximum  sensitivity 
and  is  retained  for  analysis.  Another  technique  introduces  a  90s  phase  shift 
between  onhogonal  polarizations  of  the  reference  beam.  When  the  object  beam  is 
interfered  with  each  of  the  reference  components,  the  two  resulting  interference 
intensity  signals  may  be  squared  and  then  added  thus  linearizing  the  interferometer 
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at  its  peak  sensitivity  (Vilkomeraon,  1976;  Peck,  19S3).  Finally,  optical  hetero¬ 
dyne  techniques  may  be  used  to  minimise  sensitivity  variations  caused  by  path- 
length  instabilities.  This  technique  will  be  reviewed  in  some  detail  in  following 
sections. 

It  is  usually  assumed  that  in  the  case  of  the  Michelson  interferometer  the  shot 
noise  of  the  photodetector  sets  the  fundamental  limit  to  the  interferometer’s  sensi¬ 
tivity  (Kwaaitaal,  1980).  This  may  not  always  be  the  case,  however.  Other  contri¬ 
butions  to  the  total  noise  of  the  interferometer  system  include  laser  noise,  en¬ 
vironmental  (mechanical)  noise,  and  electronic  noise.  Kwaaitaal ’s  work  indicates 
that  laser  noise  can  be  made  insignificant  if  a  single  mode  laser  is  used  and  the  fre¬ 
quency  content  of  the  detected  signal  is  above  10kHz  (Kwaaitaal,  1980).  Electron¬ 
ic  and  environmental  noise,  however,  are  not  so  easily  neglected. 


Consider  first  the  noise  associated  with  a  photodetector/amplifier  combination 
required  to  obtain  a  voltage  signal  analog  to  the  light  intensity  incident  on  the 
detector.  Two  general  classes  of  photodetectors  will  be  considered  -  photoemissive 
and  solid  state.  Specific  representatives  of  these  two  classes  are  the  photomulti¬ 
plier  tube  and  the  junction  photodiode,  respectively.  Tn  both  cases  light  incident 
on  the  active  regions  of  these  devices  (i.e.  photocathode  or  depletion  region)  li¬ 
berates  electrons  or  generates  electron- hole  pairs  which  contribute  to  current  flow 
through  the  device.  This  signal  current  may  be  described  as  follows: 


h 


(5) 


if  —  quantum  efficiency  in 


carriers/photon 
q  —  election  charge 
P  —  incident  light  intensity 
h  —  Plank’s  constant 
v  —  light  frequency. 

Owing  to  the  quantum  nature  of  light  there  is  an  associated  shot  noise  with  a  spec¬ 
tral  density  function 


S(*)"2ffq2,  (6) 

where  N  represents  the  avenge  carrier  emission  rate  (Yaiiv,  1976).  Since  the 
avenge  current/, is  oN,  the  shot  noise  power  for  a  given  frequency  bandwidth,  H, 
can  be  written: 

ifr—lqlB.  (7) 

Thus  if  shot  noise  were  the  only  source  of  noise  in  the  detector/amplifier  system, 
then  the  signal  to  noise  ratio  could  be  computed  from  Equations  5  and  7  as  fol¬ 
lows: 


SIN' 


!i 


{iiqP/hv) 

2 q(jfqP!hv)B 


The  quantum  efficiency  of  common  photocathode  material*  is  typically  about  10% 
versus  about  30%  for  a  silicon  photodiode  (  Anderson  and  McMurtry,  1966).  Thus 
if  one  considers  that  the  performance  of  a  detector/amplifier  system  is  shot  noise 
limited,  then  better  signal  to  noise  ratio  will  be  obtained  using  solid  state  photo¬ 
diode  detectors. 


Tn  a  shot  noise  limited  system  the  minimum  detectable  surface  displacement 
may  be  computed  for  the  case  of  the  Michelson  interferometer.  Equation  4  pro¬ 
vides  a  relation  for  the  power  incident  on  a  photodetector  for  a  mirror  displace¬ 
ment,  D. 


P—PdLl+Mcos  ~D), 


where  M  (modulation  factor)  —  2 


PsubO  -  a02+ar2- 


In  the  case  where  an  acoustic  wave  is  launched  on  the  mirror  surface,  the  displace¬ 
ment,  D  is  the  sum  of  a  static  component  and  a  dynamic  component.  Using 
Kwaaitaal’s  notation 


D—X+tc 


where  X  “  mirror  position 


x  — i  sin  2rrft. 


Thus  the  pcwer  incident  on  the  photodetector  is 


For  other  than  ideal  systems,  M<  1  causing  the  maximum  signal  to  noise  to  occur 
for  .if =3 A /1 6.  This  result  is  shown  graphically  in  Figure  2a.  Note  that ,  somewhat 
surprisingly,  the  best  signal  to  noise  ratio  is  not  obtained  when  x— A/8,  i-e.  where 
the  slope  of  the  signal  curve  is  greatest  (Figure  2b).  Instead  the  best  SNR  is  ob¬ 
tained  for  Af^A/lb.  One  can  reconcile  this  result  by  observing  that  the  3A/16 
point  represents  a  compromise  between  maximum  signal  gain  and  minimum  shot 
noise.  The  background  light  level  at  the  point  of  maximum  signal  slope  (X— A/8) 
is  greater  than  at  the  computed  SNR  maximum.  Therefore  the  average  photodetec¬ 
tor  current  and  consequently  the  shot  noise  will  be  greater  at  this  point.  On  the 
other  hand,  when  X— A/4,  the  shot  noise  is  at  a  minimum.  Unfortunately  the  sig¬ 
nal  gain  is  at  a  minimum  as  well.  Therefore  the  optimum  value  for  X  must  lie 
somewhere  between  A/8  and  A/4.  Using  the  argon  ion  laser  line  at  514.5nm,  and 
a  0.5mm  diameter  silicon  detector,  and  0.05mW/cm  2  output  field  intensity  typical 
of  a  holographic  interferometry  system,  the  minimum  dynamic  displacement  sensi¬ 
tivity  for  a  SNR— 1  can  be  computed  using  Equation  10.  For  X— 3A/16,  M— 0.9, 
and  a  1Hz  bandwidth, 

x— 4.71xl(Tflm. 

Note  that  this  is  about  500a  and  is,  therefore,  far  less  sensitive  than  desirable. 
Fonunately,  sensitivity  may  be  improved  by  increasing  laser  illumination,  holo¬ 
gram  efficiency,  detector  area,  and  by  using  heterodyne  detection  all  of  which  tend 
to  increase  the  light  intensity  incident  on  the  photodetector.  Rear  in  mind  that  this 
rather  disappointingly  poor  sensitivity  is  largely  the  result  of  the  low  light  intensity 
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Figure  2:  a)  Signal  to  Noise  Ratio  and  b)  Photocurrent 
as  a  function  of  path length  difference,  X,  in 
a  Kichelson  interferometer. 


figure  which  was  selected  to  represent  actual  intensities  encountered  in  holographic 
interferometry. 


Heterodyne  Interferometry 


Heterodyne  interferometry  provides  another  alternative  to  path  stabilization  for 
immunity  to  environmental  noise  sources.  In  addition,  the  application  of  hetero¬ 
dyne  techniques  to  holographic  interferometry  allows  for  interpolation  between 
fringes  to  obtain  sensitivity  to  static  displacements  comparable  to  real  time  dynam¬ 
ic  detection  sensitivities  obtained  by  classical  interferometers.  The  basic  difference 
between  heterodyne  and  homodyne  interferometry  is  that  the  two  wave  fronts  to 
be  interfered  in  the  output  plane  of  the  interferometer  are  shifted  in  frequency  re¬ 
lative  to  each  other  for  heterodyne  and  not  for  homodyne.  The  light  intensity  in¬ 
cident  on  the  photodetector  may  be  computed  in  a  manner  similar  to  that  used  in 
the  development  of  Equation  3. 


Then  the  intensity,  I,  at  the  detector  is: 


Hu  x+u  JMi/V +u\ xu  2) 


•A  f  -K4  }  ~h2A  jA^cosKoj^ — 


(11) 


The  photodetector  signal,  then,  is  an  offset  sinusoid  whose  intensity  varies  at  the 
difference  frequency,  /— 2Tr(<u1— **2).  When  used  in  a  Michelson  interferometer 
configuration,  a  dynamic  displacement  of  the  object  mirror  results  in  a  phase  shift 
in  the  detected  signal.  Therefore,  displacements  are  detected  using  electronic 
phase  measurement  techniques  which  can  be  made  insensitive  to  signal  amplitude 
variations. 

Several  techniques  have  been  used  to  detect  and  process  the  phase  shifted  pho¬ 
todetector  signal.  Frequency  demodulation  has  been  used  to  measure  the 
Doppler-shifted  output  signal  (Whitman, .  Laub,  and  Hates,  1968).  The  doppler 
shifted  frequency  is  proportional  to  the  velocity  of  the  surface  motion  and  there¬ 
fore  electronic  integration  must  be  performed  in  order  to  obtain  displacement  in¬ 
formation.  The  difficulties  in  accurately  performing  the  required  integration  partic¬ 
ularly  at  frequencies  very  near  the  difference  frequency  make  Doppler  frequency 
detection  a  less  attractive  option  for  heterodyne  signal  processing  (Palmer  and 
Green,  1979).  FM  detection  using  a  phase  locked  loop  has  been  used  to  success¬ 
fully  detect  high  frequency  vibration  and  pulsed  displacements  (Rudd,  1980).  In 
such  a  system  an  electronically  phase  locked  loop  is  locked  to  the  heterodyne  sig¬ 
nal  from  the  photodetector.  Phase  variations  occurring  at  a  rate  below  the  loop 
response  rate  (typically  100kHz)  are  tracked  by  the  loop.  Thus  the  output  of  the 
phase  detector  remains  constant  for  low  frequency  phase  changes.  Since  the  fre¬ 
quency  content  of  most  environmental  noise  sources  is  well  below  100kHz,  the 
loop  provides  a  level  of  noise  immunity  similar  to  that  obtained  with  path  length 


stabilized  interferometers.  For  vibrational  disturbances  above  100kHz,  the  error 
signal  present  at  the  output  of  the  phase  detector  provides  a  voltage  signal  propor¬ 
tional  to  the  velocity  of  the  vibrational  disturbance.  Another  scheme  also  sensitive 
to  the  velocity  of  the  test  object’s  surface  displacement  uses  coherent  detection  of 
the  phase  shifted  signal  (DeLaRue,  1972).  This  technique  works  particularly  well 
when  the  frequency  of  the  acoustic  exciting  source  is  well  known.  In  this  case  the 
signal  to  be  detected  is  shifted  to  a  sideband  of  the  detector  carrier  frequency.  By 
using  the  earner  signal  as  a  local  oscillator,  immunity  to  environmentally  induced 
path  length  instability  is  obtained.  This  optically  derived  local  oscillator  is  then 
mixed  with  one  of  the  signal  sidebands  to  produce  the  desired  output  signal.  In 
another  coherent  detection  scheme  for  phase  demodulation,  the  local  oscillator  is 
simply  the  input  signal  delayed  by  a  known  amount,  typically  =  500ns  (Jacobson, 
1982).  By  delaying  the  local  oscillator  by  some  integral  number  of  wavelengths 
plus  A/4,  it  can  be  made  to  be  in  quadrature  with  the  heterodyne  signal.  Coherent 
detection  using  these  two  signals  and  subsequent  filtering  can  be  used  to  to  provide 
a  direct  measure  of  phase.  Again  low  frequency  noise  due  to  pathiength  instability 
may  be  filtered  out. 

When  changes  in  the  detected  phase  difference  occur  slowly,  direct  phase 
measurement  of  the  detected  carrier  frequency  may  be  used.  To  do  so  requires  a 
local  oscillator  source  coherent  with  the  detector  signal.  If  the  relative  phase 
difference  between  the  two  is  measured,  then  any  change  in  the  optical  path 
difference  between  the  reference  and  object  beams  will  cause  a  corresponding 
change  in  the  measured  phase  difference.  Electronic  phase  difference  measure¬ 
ments  are  quite  easily  and  most  commonly  made  using  electronic  instrumentation 


which  produces  an  analog  output  voltage  proportional  to  the  time  between 
corresponding  zero  crossings  of  the  two  input  signals  expressed  as  a  fraction  of  the 
signal  period.  This  type  of  phase  meter  is  capable  of  providing  accuracies  better 
than  0.1°  (Mastner,  1980)  thus  giving  a  potential  optical  path  difference  sensitivity 
of  (A/2X0.1/360)  —  X  x  1.4  x  10-4  or  about  0.07nm.  For  the  work  reported  here, 
the  relative  phase  differences  between  the  signal  and  reference  oscillator  were 
derived  from  the  outputs  of  a  quadrature  detection  lock-in  amplifier.  Using  a 
lock-in  amplifier  also  reduces  the  effects  of  variations  in  input  signal  amplitude  en¬ 
countered  with  other  phase  measurement  schemes  (Ineichen,  1977). 

Coherent  quadrature  detection  for  relative  phase  measurements  is  shown 
schematically  in  Figure  3.  The  input  signal  is  mixed  both  with  the  local  oscillator 
and  with  a  90s  phase  shifted  version  of  the  local  oscillator. 

Signal,  S  “  Acos(«ut  4-  <&)  +  N 

Local  oscillator,  R  (  —  Hcosont 
Shifted  L.O.,  R  2  —  Hsinwt. 

Output  from  first  mixer  is 

O  ]-mSR  \—ABcos  toteos  ( wt  -Hi  y+NBcas  aW 

—^-)ccs<b+NBcostU. 


Output  from  the  second  mixer  is 


O  f-SR  f“ABsin  <wt  cosiati  -H6  'H-NBsin  tut 
— (  )sin  (  2arf  -Hi  )—<  )stn  <jr±NBan  td. 

Note  that  by  low  pass  filtering  the  outputs  O  (  and  O  only  the  second  term  in 
each  expression  is  retained  (assuming  that  A  tB  and  4  are  slowly  varying  func¬ 
tions  of  time).  Note  also  that  the  additive  noise,  N,  originally  present  in  the  input 
signal  is  eliminated  by  this  filtering  process  as  wed.  Tn  fact  it  is  by  this  detection 
and  filtering  process  that  lock-in  amplifiers  may  detect  and  amplify  signals  in  noise 
backgrounds  as  much  as  1000  times  greater  than  the  signal  amplitude!  The  signals 
remaining  after  the  filtering  process  then  are  as  follows: 


OFi^~  sin* 

By  taking  the  arctan  of  (OF -J0F\),  the  relative  phase  difference,  </»,  is  computed. 
Some  lock-in  amplifiers  are  available  which  compute  arctaniOF JOF ,)  by  an  ana¬ 
log  process  in  real  time.  Again  owing  to  the  coherent  detection  and  filtering  pro¬ 
cess,  the  ratio  (OF-JOF t)  and  0,  therefore,  is  independent  of  the  signal  amplitude, 
A. 

Judicious  choice  of  the  local  oscillator  source  can  reduce  the  effects  of  environ¬ 
mental  noise  on  the  accuracy  of  the  relative  phase  measurement.  The  most  com¬ 
monly  used  technique  for  providing  optical  frequency  shifting  in  heterodyne  inter¬ 
ferometry  is  the  use  of  acousto-optic  modulators,  or  Bragg  cells.  Other  techniques 
which  have  been  used  include  rotating  diffraction  gratings,  rotating  wave  plates. 


and  electro-optical  "serrodyning"  (Hard,  1978;  Aleksoff,  1971).  Using  a  Bragg  cell, 
the  first  order  diffracted  beams  exiting  the  cell  are  shifted  into  upper  and  lower 
sidebands  around  the  optical  frequency  at.  Light  in  either  of  the  diffracted  orders 
varies  in  frequency  from  the  undiffracted  beam  by  ±fD,  where  f D  is  the  RF  driv¬ 
ing  frequency  of  the  cell.  The  frequency  shift  is  the  result  of  a  phenomenon 
kmvn  as  the  Debye-Sears  effect  (Adler,  1967).  Typically  the  RF  drive  frequency 
for  commercial  Bragg  cells  is  40  or  80MHz.  Therefore  if  an  80MHz  Bragg  cell  is 
used  to  perform  the  frequency  shifiting  in  a  heterodyne  interferometer,  the 
difference  frequency  signal  from  the  photodetector  will  be  equal  to  the  80 MHz 
drive  frequency.  While  a  high  frequency  carrier  may  be  desirable  in  real-time 
detection  of  acoustic  displacements,  it  is  not  necessary  for  the  detection  of  slowly 
varying  changes  in  optical  path  length.  Indeed  the  design  and  construction  of 
amplifiers  and  *«nrim«ri  processing  electronics  becomes  somewhat  more  compli¬ 
cated  in  the  higher  frequency  regimes. 

A  more  manageable  carrier  frequency  is  obtained  by  passing  the  laser  beam 
through  two  Bragg  cell  modulators  driven  at  slightly  different  frequencies  as  shown 
in  Figure  4.  Light  diffracted  into  the  + 1  order  of  the  first  ceil  is  frequency  shifted 
upward  by  80MHz.  This  first  order  diffracted  beam  passes  subsequently  through 
the  second  Bragg  cell.  The  -1  order  beam  exiting  the  second  cell  is  now  shifted 
down  in  frequency  by  80.1MHz.  Therefore  the  net  frequency  shift  caused  by  pass¬ 
ing  the  beam  through  both  cells  is  only  0.1MHz  (100kHz).  The  resulting  100 kHz 
difference  frequency  detected  at  the  output  of  the  heterodyne  interferometer  is  far 
easier  to  amplify  and  process  than  would  be  the  80MHz.  The  diffraction  efficiency 
of  commercially  available  Bragg  cell  modulators  using  lithium  niobate  is  better 


Figure  4:  Obtaining  a  100kHz  frequency- 
shift  using  two  Eragg  cells 


than  80%.  Consequently,  the  overall  efficiency  of  this  frequency  shifting  process  is 
on  the  older  of  65%.  Unfortunately  it  is  impractical  to  use  just  a  single  cell  at 
100k Hz  driving  frequency  because  the  diffraction  efficiency  would  be  quite  poor, 
and  more  importantly,  an  extremely  small  diffraction  angle  of  the  100kHz  shifted 
components  makes  it  difficult  to  spatially  resolve  the  diffracted  and  undiffracted 
orders.  There  is,  however,  one  very  convenient  feature  of  practical  importance  as¬ 
sociated  with  using  two  Bragg  cells  for  frequency  shifting.  That  is  the  fact  that  both 
cells  may  be  driven  at  the  same  shift  frequency  resulting  in  no  net  frequency  shift 
nor  causing  any  appreciable  change  in  the  optical  path  or  angle  through  'he  cells. 
Therefore  the  heterodyne  interferometer  may  be  used  in  a  homodyne  mode 
without  the  need  for  realiignment  of  the  beam  in  the  shifted  leg.  This  feature  is 
particularly  useful  in  initial  allignment  of  the  heterodyne  interferometer. 

Most  of  the  discussion  to  this  point  has  been  with  regard  to  the  use  of  inter¬ 
ferometers  to  detect  and  measure  dynamic  displacements  at  a  single  point  on  the 
surface  of  an  object  minor.  In  fact,  the  research  reported  in  this  paper  also  ad¬ 
dresses  the  problem  of  detecting  and  mapping  dynamic  surface  displacements  asso¬ 
ciated  with  traveling  surface  acoustic  waves.  However,  unlike  the  inierferometric 
techniques  which  have  been  discussed  to  this  point,  heterodyne  holographic  inter¬ 
ferometry  is  not  limited  to  real  time  surface  displacement  measurements  at  a  single 
point  on  the  test  object.  Instead,  the  surface  displacements  over  the  entire  field  of 
view  of  the  test  object  are  recorded  at  a  single  moment  in  time  using  pulsed  holog¬ 
raphy.  Subsequently,  heterodyne  interferometric  techniques  are  applied  to  the 
reconstructed  holographic  image  to  provide  details  of  the  surface  deformation  over 
the  entire  field  of  view.  This  later  step  is  quite  similar  to  broad  beam  heterodyne 


interferometry  which  has  been  used  to  measure  surface  flatness  of  optical  com¬ 
ponents,  for  example. 

In  broad  beam  interferometry,  the  illuminating  beams  are  expanded  and  col¬ 
limated  so  that  large  surface  areas  may  be  interferometrically  compared.  If  the  ob¬ 
ject  mirror  is  tilted  or  deformed  relative  to  the  reference  mirror,  a  fringe  pattern  is 
generated  in  the  output  plane  as  shown  in  Figure  5.  The  fringes  map  out  contours 
of  constant  displacement  with  a  contour  interval  of  X/2.  When  a  relative  frequen¬ 
cy  shift  is  introduced  between  the  reference  and  object  beams,  a  detector  placed  in 
the  output  field  generates  a  sinusoidal  signal  at  the  shift  frequency.  Direct  meas¬ 
urement  of  the  phase  difference  between  this  heterodyne  signal  and  a  local  refer¬ 
ence  oscillator  may  then  be  used  to  interpolate  between  fringes  providing  a  direct 
measure  of  out  of  plane  displacement  as  a  function  of  position  on  the  object  mirror 
surface.  By  scanning  the  detector  over  the  output  field,  a  map  of  surface  topogra¬ 
phy  may  be  generated.  Such  real-time  broad  beam,  heterodyne  interferometry  has 
been  used  to  measure  flatness  and  curvature  in  mirrors,  silicon  wafers,  and  optical 
components  (Massie,  1979). 

Broad  beam  interferometry  provides  yet  another  possibility  for  reducing  noise 
resulting  from  path  instabilities  in  interferometry.  A  reference  detector  may  be 
placed  at  any  fixed  point  in  the  output  field.  The  sinusoidal  signal  from  this  refer¬ 
ence  detector  may  be  used  as  a  local  oscillator  signal  against  which  the  phase  of  the 
signal  from  the  scanning  detector  may  be  measured.  Low,  frequency  vibrational 
disturbances  which  may  affect  the  optical  pathlengths  of  either  the  reference  or  ob¬ 
ject  leg  of  the  interferometer  will  cause  phase  shifts  in  both  the  reference  and 


scanning  detector  signals.  Since  both  detectors  are  placed  in  the  same  output 
plane,  the  phase  shifts  resulting  from  vibration  will  be  equal  in  both  signals. 
Therefore  the  easured  phase  difference  between  the  two  signals  will  be 
unaffected  by  the  vibrational  disturbance. 

Computing  the  noise  limited  sensitivity  of  a  broad  beam  heterodyne  inter¬ 
ferometer  requires  an  understanding  of  the  effect  of  signal  to  noise  ratio  (SNR)  on 
the  performance  of  the  system’s  phase  meter.  Assuming  that  the  detector 
bandwidths  are  less  than  the  heterodyne  carrier  frequency,  then  a  phase  error,  A 
using  a  zero-crossing  phase  meter  is 


where  r—  integration  time 
f  —  heterodyne  frequency 
(Ineichen,  1977). 

Tf  phase  fluctuations  from  vibration  and  other  environmental  effects  occur  at  fre¬ 
quencies  outside  of  the  system  bandwidth,  then  the  SNR  can  be  assumed  to 
depend  only  upon  shot  noise. 

To  compute  the  shot  noise  limited  SNR,  recall  from  Equation  1 1  that  the  light 
intensity  incident  on  the  photodetector  is 


for  M— 1A  tA  -JA  2  +v422 
and 


The  avenge  power,  P<>  contributes  to  the  shot  noise  current  while  the  signal 
current  varies  with  the  cosine  term: 
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The  phase  error  resulting  from  a  shot  noise  limited  system  which- incorporates  a 
zero-crossing  detector  phase  meter  is 


A* 


_ 1 _ 

yJ(i)M1P(J^h  vB)fr 


(12) 


Using  the  same  values  for  tj,  M,  P0>  a™!  B  as  used  earlier  for  the  example  calcula¬ 
tion  for  the  path-stabilized  interferometer  (i.e.  rj  ”  0.3,  M  ”  0.9,  Ptf-lO-7  watts, 
and  B  —  1Hz)  and  using  an  Ar+  laser  (X  —  05145micron)  and  a  zero-crossing 


phase  meter  where  / r— 1,  the  phase  error,  A6,  for  a  heterodyne  interferometer  is 
8x  1CT6  radians  corresponding  to  a  displacement  sensitivity  of  about  0.003a.  Note 
that  this  is  some  5  orders  of  magnitude  better  than  was  obtained  for  the  path  sta¬ 
bilized  interferometer. 


Holographic  Interferometry 


To  this  point,  all  discussion  of  interferometry,  path  stabilization,  and  detection 
has  been  illustrated  by  consideration  of  the  Michelson  interferometer.  Holographic 
interferometry  incorporates  many  of  the  same  operating  principles  but  reduces 
some  restrictions  imposed  by  classical  interferometry,  particularly  regarding  the 
types  of  specimens  which  may  be  examined.  For  any  form  of  interferometry  to  be 
of  practical  use,  the  interference  fringe  pattern  in  the  output  plane  must  be  of  rela¬ 
tively  low  order  so  that  individual  fringes  may  be  resolved.  Obviously  if  the  fringe 
rate  is  too  high,  there  will  be  no  fringe  paitern  observed  visually  or  detected  with  a 
photodetector  of  finite  area  with  dimensions  larger  than  one  fringe  period.  This 
necessity  for  a  low  fringe  rate  imposes  the  restriction  that  the  two  wavefronts  to  be 
interfered  must  be  very  nearly  identical.  Tn  the  case  of  a  classical  interferometer, 
nominally  identical  wavefronts  are  reflected  from  two  flat  mirrors  with  one  mirror 
serving  as  a  reference.  When  the  object  mirror  is  moved,  tilted,  or  deformed, 
differences  between  the  reference  and  object  wavefronts  occur  which  manifest 
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themselves  in  the  observed  interference  pattern.  If  the  object  mirror  is  replaced  by 
a  diffusely  reflecting  surface,  a  very  complicated  and  high  fringe  rate  pattern  ap¬ 
pears  in  the  output  plane.  The  source  of  this  complex  pattern  is  clear  if  one  con¬ 
siders  that  a  diffusely  reflecting  surface  is  made  up  of  a  mosaic  of  small  mirror-flat 
portions  at  random  angles  each  reflecting  a  portion  of  the  illuminating  beam  back 
to  the  output  plane.  Therefore  the  fringe  pattern  in  the  output  plane  is  the  result 
of  interference  of  the  light  reflected  from  each  of  the  individual  reflectors  with  the 
reference  beam  as  well  as  with  each  other.  (It  is  the  interference  pattern  from  the 
interaction  of  the  light  from  the  individual. surface  reflectors  with  each  other  which 
is  called  speckle.  The  effects  of  speckle  on  the  accuracy  of  holographic  inter¬ 
ferometry  measurements  will  be  discussed  in  some  defail  shortly.)  By  taking  ad¬ 
vantage  of  the  "mosaic"  nature  of  diffuse  surfaces,  one  may  still  use  classical  inter¬ 
ferometry  to  transduce  surface  displacements.  This  is  often  done  by  focusing  the 
object  illuminating  beam  down  to  a  small  point  such  that  at  the  point  of  illumina¬ 
tion  the  local  variation  of  the  object  surface  is  negligible.  The  light  is  reflected 
directly  back  from  the  test  point  and  recollimated  by  the  focusing  objective  lens. 
Subsequent  interference  of  this  object  beam  with  the  reference  takes  place  as 
though  the  beam  were  reflected  from  a  flat  mirror. 

There  are  a  number  of  applications,  however,  where  it  is  desirable  to  inter- 
feroraetrically  examine  an  entire  surface.  As  discussed  previously,  broad  beam  in¬ 
terferometry  may  be  used  for  objects  with  specularly  reflecting  surfaces  which  are 
either  flat  or  of  a  known  simple  geometry.  In  such  a  system,  a  "master"  object 
against  which  a  test  object  is  to  be  compared  is  used  in  place  of  a  reference  mirror. 
Usually  wavefront  correcting  lenses  are  placed  in  each  leg  of  the  interferometer  so 
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that  the  extent  of  the  beams  reflecting  from  each  object  can  be  controlled  to  stay 
within  the  physical  dimensions  of  the  components  comprising  the  interferometer. 
The  use  of  conventional  spherical  lenses  restricts  the  nominal  shapes  of  the  test 
object  surfaces  to  be  either  flat,  convex,  or  concave. 

Interferometric  analysis  of  more  complex,  diffusely-reflecting  surfaces  requires 
that  a  reference  wavefront  be  generated  which  matches  that  of  the  test  object  sur¬ 
face  One  might  consider  the  possibility  of  manufacturing  a  roaster  object  to  re¬ 
place  the  reference  reflector  in  a  Michelson  interferometer  as  described  above.  It 
is  important  to  remember,  however,  that  in  order  for  such  an  approach  to  work, 
every  dimension  of  the  master  object  including  each  minute  reflector  in  the  "mosa¬ 
ic'*  comprising  the  object  surface  must  match  that  of  the  test  object  to  be  exam¬ 
ined.  It  is  certainly  impractical,  if  not  impossible,  to  produce  reference  surfaces 
with  the  required  precision. 

Holography  provides  one  solution  to  the  problem  of  generating  a  master  wave- 
front  for  interferometric  comparison  of  objects  with  arbitrary  surface  shape  and 
texture.  While  some  work  has  been  done  adapting  holography  for  interferometric 
comparison  of  different,  nominally  identical  objects  (Neumann,  1981),  by  far  the 
most  common  applications  of  holography  in  interferometry  is  in  real-time  or 
double-exposure  interferometry  where  the  object  surface  is  compared  at  two 
different  times.  By  interfering  the  wavefronts  which  reflect  from  the  object  both  in 
a  stressed  and  an  unstressed  state,  for  example,  one  may  observe  deformations 
which  occur  in  the  object  as  a  result  of  applied  stress.  If  the  image  of  an  object  as 
recorded  on  a  hologram  is  interfered  with  the  image  of  the  object  at  a  slightly 


different  longitudinal  magnification,  the  object  appears  superimposed  with  interfer¬ 
ence  fringes  forming  contours  of  constant  surface  elevation. 

The  operating  principles  of  holography  may  be  explained  using  either  wave- 
front  interference  theory  front  geometrical  optics  or  alternatively  using  scalar 
diffraction  theory.  It  is  the  latter  approach  which  will  be  pursued  here.  In  scalar 
diffraction  theory,  often  called  Fourier  Optics  (Goodman,  1968),  an  optical  wave- 
front  may  be  described  as  a  complex  quantity  whose  magnitude  and  phase  varies  as 
a  function  of  x  and  y  in  a  plane  normal  to  the  optical  axix,  z.  A  phasor  notation  of 
sorts  may  be  used  to  explicitly  specify  the  propagation  direction  of  a  wavefront  as 
the  angle  between  the  propagation  direction  and  the  optical  (z)  axis  as  follows 
(Collier,  1971). 

(5<Oejr(t'>~OeJ* 

where  <J*-/(0)a0. 

Looking  at  the  optical  set  up  shown  in  Figure  6,  note  that  two  wavefronts  intersect 
at  the  film  plate.  One  wavefront  results  from  the  reflection  of  the  illuminating 
beam  by  the  object  surface  and  will  be  referred  to  as  the  object  wavefront,  (5 .  The 
second  beam  is  unmodified  by  reflection  and  is  called  the  reference  beam,  R. 
Defining  the  optical  axis  to  be  perpendicular  to  the  film  plate: 

R~ReJ* 

6  -Oe°~0. 


Thus  the  optical  field  falling  on  the  film  plate,  A  —Rr1  *+0 .  Assuming  that  the 
transmittance,  T,  of  the  exposed  and  subsequently  developed  film  is  proportional 
to  the  field  intensity,  I,  then: 

HR+OKR'+O') 

HReJ  *40  X/fe_>*40* ) 

^R\^Oy-OR'  e~j*+0'  ReJ*.  (13) 

Upon  reillumination  by  the  reference  beam  as  shown  in  Figure  7,  the  hologram 
diffracts  the  incident  light  into  three  beams  corresponding  to  the  three  bracketed 
terms  below. 

(14) 

Analyzing  Equation  14,  the  first  term  describes  the  intensity  variations  made  to  the 
reconstructing  reference  beam  as  it  passes  through  the  hologram  plate  and  contin¬ 
ues  to  propagate  at  an  angle,  0,  to  the  optical  axis.  The  third  term  contains  the 
conjugate  image  of  the  object  and  leaves  the  film  plate  at  an  angle  20  relative  to 
the  optical  axis.  It  is  the  second  term  of  Equation  14  that  is  of  importance  in  holo¬ 
graphic  interferometry-  Note  that  this  term  describes  a  wavefront  which  pro¬ 
pagates  along  the  optical  axis  (4>— 0-0)  and  is  the  product  of  the  original  object 
wavefront  multiplied  by  the  squared  magnitude  of  the  reference  wavefront.  In  the 
case  being  considered  here,  the  reference  wave  is  a  uniform,  collimated  beam  such 
that  its  squared  magnitude  is  simply  a  scalar  constant  independent  of  x  and  y.  Thus 
a  reconstructed  replica  of  the  original  object  wavefront  propagates  along  the  optical 
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Figure  7:  Wavefront  reconstruction  by  holography. 
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•xu  such  that  an  observer  portioned  as  shewn  in  Figure  7  would  see  a  virtual  im¬ 
age  of  the  otgect  appearing  in  its  original  location  and  in  3-D. 

The  holographs caiiy  recorded  image  of  the  object  may  new  be  used  as  a  'mas¬ 
ter'  in  an  interferometer  as  shewn  in  Figure  8.  The  recording  and  playback 
geometries  shown  in  Figures  6,  7,  and  8  were  chosen  to  illustrate  the  similarity  of 
holographic  interferometry  to  the  classical  Michelson  interferometer.  In  practice, 
however,  off-axis  illumination  of  the  test  object  is  used  and  the  hologram  serves  as 
both  a  beam  splitter  and  reference  image  source.  Off-axis  illumination  does,  how¬ 
ever,  complicate  fringe  analysis  somewhat.  Remember  that  interferometric  fringe 
formation  is  the  result  of  optical  pathlength  differences  which  occur  between  the 
reference  and  object  wavefronts.  For  coaxial  illumination,  i.e.  when  the  illumina¬ 
tion  and  viewing  axes  are  the  same,  one  fringe  corresponds  to  an  object  displace¬ 
ment  of  A/2.  For  off-axis  illumination,  however,  the  sensitivity  to  object  displace¬ 
ment  is  somewhat  reduced.  For  analytical  purposes  a  sensitivity  vector,  K,  is 
defined  in  the  direction  of  the  bisector  of  the  illuminating  and  viewing  axes  (see 
Figure  9)  (Vest,  1979).  The  sensitivity  vector  has  a  magnitude: 

(X|-2Aroj0 
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Figure  8:  Holographic  interferometry  using 
coaxial  illumination. 
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Figure  9:  Sensitivity  vector,  K,  for  off-axis 
holographic  interferometry. 


For  an  object  displacement  D,  the  observed  intensity  in  the  output  plane  of  the 
holographic  interferometer  varies  in  a  manner  similar  to  that  described  in  Equation 
4  as  follows: 

I~2A\  1-hcaj  K  D)  (15) 

—2A  ^  14-C05  ( IkDcos  0cos<6 )  ] , 

where  is  the  angle  between  the  sensitivity  vector  and  the  displacement  direction, 
D.  Another  advantage  of  holography  is  now  apparent.  Using  Equation  IS  the 
component  of  D  in  the  K  direction  may  be  determined  for  any  point  on  the  object 
surface.  By  selecting  three  different  viewing  angles,  the  projection  of  D  on  three 
different  K  vectors  may  be  computed  and  subsequently  processed  to  completely 
specify  D  in  magnitude  and  direction. 

Double  exposure  holography  eliminates  the  need  to  precisely  replace  and  align 
the  hologram  in  the  optical  set-up  after  its  exposure  and  development.  Tn  this  case, 
both  the  reference  and  object  images  are  reconstructed  by  the  hologram.  The  opti¬ 
cal  arrangement  for  recording  a  double  exposure  hologram  can  be  identical  to  that 
used  in  real  time  holographic  interferometry.  Two  exposures  of  the  object  image 
are  recorded  with  a  change  either  in  magnification  or  applied  stress  being  made 
between  exposures.  The  alignment  of  the  reference  and  object  images  is  preserved 
in  the  hologram  so  that  upon  reconstruction  with  a  single  reference  beam,  both 
images  are  superimposed  and  therefore  interfere.  One  variation  of  this  double¬ 
exposure  geometry  is  important  in  heterodyne  applications.  This  variation  is 
shown  in  Figure  10  where  two  spatially  separated  reference  beams  are  used  to 
record  and  reconstruct  the  holographic  images.  One  reference  beam  and  then  the 
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Figure  10:  Dual  reference  beam  holography 
for  heterodyne  interferometry. 


other  is  used  in  making  the  two  exposures  respectively.  On  playback,  both  beams 
are  used  to  simultaneously  reconstruct  the  reference  and  object  images.  If  one  of 
the  reconstrudng  reference  beams  is  shifted  in  frequency  relative  to  the  other, 
then  heterodyne  analysis  techniques  similar  to  those  used  in  bread  beam  hetero¬ 
dyne  interferometry  may  be  used. 


Heterodyne  Holographic  Interferometry 


Heterodyne  analysis  of  dual  reference  beam,  double  exposure  holograms 
recorded  as  described  above,  permits  interpolation  between  the  normally  observed 
interferometric  fringes.  By  interpolating  between  fringes  both  the  sensitivity  and 
dynamic  range  of  holographic  interferometry  may  be  greatly  improved-  When  first 
introduced  in  1973,  interpolation  to  better  than  1/100  of  a  fringe  was  demonstrat¬ 
ed  (Dandliker,  1973).  In  1980,  Dandliker  published  a  rather  complete  review  pa¬ 
per  on  heterodyne  holographic  interferometry  which  predicts  practical  sensitivity 
limits  of  better  than  1/1000  of  a  fringe  (Dandliker,  1980).  Indeed  Dandliker  and 
his  associates  dominate  the  literature  from  1973  to  1980  on  this  topic.  Consequent¬ 
ly  much  of  the  theoretical  development  to  be  presented  here  is  discussed  in  greater 
detail  in  his  1980  publication. 


In  many  respects,  the  principles  of  heterodyne  holographic  interferometry 


parallel  those  of  broad,  beam  heterodyne  interferometry  discussed  earlier.  In  the 
holographic  case,  however,  the  two  wavefronts  to  be  interfered  are  reconstructed 
from  the  hologram  plate  rather  than  being  reflected  from  reference  and  object  mir¬ 
ror  surfaces.  By  constructing  the  hologram  with  two  reference  beams  as  shown  in 
Figure  10,  the  two  reconstructed  image  wavefronts  interfere  only  when  both  refer¬ 
ence  beams  illuminate  the  developed  hologram  simultaneously.  When  one  refer¬ 
ence  beam  is  shifted  in  frequency  relative  to  the  other,  the  two  reconstructed  im¬ 
age  wavefronts  are  likewise  shifted  in  frequency  causing  a  corresponding  variation 
in  the  image  intensity.  Mathematically,  the  transmittance  of  the  dual  reference 
beam,  double  exposure  hologram  can  be  expressed  as  the  superposition  of  the 
transmittances  which  would  result  from  each  exposure  independently  as  shown 
previously  in  Equation  13. 

T a  (16) 

where  the  subscripts,  1  and  2,  refer  to  the  first  and  second  exposures,  respectively. 
Again,  only  the  last  two  terms  will  give  rise  to  wavefronts  which  will  propagate 
along  the  optical  axis  to  the  output  plane. 

To  perform  heterodyne  analysis  for  fringe  interpolation,  the  two  reference 
beams  used  to  reconstruct  the  object  wavefront,  Ot  and  02,  must  be  at  slightly 
different  frequencies.  Let 
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When  together  they  illuminate  the  hologram  whose  transmittance  is  as  shewn  in 
Equation  16,  the  light  intensity  failing  on  the  output  plane  of  the  interferometer  is 
obtained  as  follows: 


If  R  ,  and  /?  2  are  plane  waves  of  unit  amplitude  then 

HO  \eJm*\2. 

Note  that  only  the  last  two  terms  of  this  expression  correspond:  to  wavefronts 
which  propagate  exactly  down  the  optical  axis  toward  the  output  plane.  The  first 
two  terms  form  images  which  lie  to  either  side  of  the  optical  axis  propagating  at  an 
angle  equal  to  that  between  the  reference  beams.  Therefore  considering  only  the 
last  two  terms  above  and  noting  that 

Ot-oieJ* 


Then 
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(17) 


This  is  exactly  the  same  form  as  obtained  in  Equation  11.  Therefore  sub-fringe  in¬ 
terpolation  between  interferometric  fringes  may  be  accomplished  at  any  point  in 
the  image  by  scanning  a  detector  over  the  image  plane  and  applying  any  of  the 
phase  detection  schemes  discussed  previously. 

Unlike  the  case  of  simple  heterodyne  interferometry,  the  resolution  of  hetero¬ 
dyne  holographic  interferometry  is  limited  both  by  shot  noise  and  speckle  noise.  If 
a  zero-crossing  phase  meter  is  used  to  measure  the  phase  difference  between  sig¬ 
nals  from  two  detectors  located- in  the  output  plane,  then  the  shot  noise  limited 
phase  sensitivity  is  the  same  as  that  of  broad  beam  heterodyne  interferometry  as 
expressed  in  Equation  12.  To  determine  the  effect  of  speckle  on  the  accuracy  of 
phase  measurements  it  is  helpful  to  understand  some  of  the  statistical  properties  of 
speckle  The  discussion  which  follows  is  not  intended  to  provide  a  complete 
description  of  speckle  phenomena  but  rather  will  attempt  to  develop  those  proper¬ 
ties  of  speckle  which  directly  affect  the  accuracy  of  heterodyne  holographic  inter¬ 
ferometry. 

Speckle  is  the  name  given  to  the  phenomenon  which  gives  a  granular  appear¬ 
ance  to  diffusely  reflecting  (or  transmitting)  objects  illuminated  by  a  coherent 
source.  The  origin  of  speckle  is  quite  easily  explained  when  one  considers  that  a 
diffusely  reflecting  (rough)  surface  can  be  considered  to  be  composed  of  a  mosaic 
of  mirror-flat  sections  whose  dimensions  are  on  the  order  of  the  wavelength  of  the 


illuminating  light  and  whose  faces  are  oriented  randomly  about  the  average  surface 
normal.  At  any  point  in  front  of  such  an  object,  light  from  the  illuminating  source 
is  scattered  by  a  number  of  surface  reflectors  creating  a  complex  interference  pat¬ 
tern.  In  general,  this  speckle  pattern  is  observed  using  an  imaging  system  such  as 
the  eye  or  a  camera  so  that  there  is  a  finite  aperture  through  which  the  light  must 
pass  before  it  is  focused  on  the  image  plane.  This  finite  aperture  limits  the  resolu¬ 
tion  of  the  imaging  system  thus  establishing  the  speckle  size  at  that  limit.  Conse¬ 
quently  for  a  square  or  a  circular  aperture,  the  speckle  size  is  Xz/d  or  1.22Xz/D, 
respectively,  where  z  is  the  distance  from  the  aperture  to  the  image  plane  and  D  is 
the  aperture  width  (Jenkins  and  White,  1950).  The  results  of  this  rather  intuitive 
argument  for  determining  speckle  size  agree  with  those  obtained  by  more  rigorous 
analysis  (Dandliker,  1980;  Vest,  1979;  Goodman,  1975).  Speckle  size,  as  will  be 
shown  shortly,  is  an  important  factor  in  determining  fringe  contrast  of  holographic 
interferograms. 

As  discussed  earlier,  the  only  practical  way  to  interfere  wavefronts  reflecting 
from  diffuse  object  surfaces  is  through  holographic  interferometry.  In  the 
mathematical  description  of  the  two  reflected  fields,  a  roughness  term,  p,  must  be 
added  to  account  for  the  microscopic  variations  in  surface  topography  which  give 
rise  to  speckle.  Therefore  the  two  wave  fields  are  expressed  as  follows: 

/i2-a2(x0)pi(xo). 


where  a ,  and  o2  represent  the  macroscopic,  or  smooth,  variations  in  the  reflected 
wavefront  and  xorfxoj'o)  positions  in  the  object  plane.  On  the  output  plane  of  a 


[inear  space- in vanent  imaging  system,  the  complex  amplitude  of  each  wavefront 
may  be  computed  by  convolution  with  the  impulse  response  of  the  imaging  sys¬ 
tem,  h(x ): 


3,(x,  >-Ja  ](xo)pi(x0)h  {Xj  —Mxo)dx0  (18) 

B^x,  >-J*i a  2(x0)p2(xo)/»  (Xj  -Mx0)dx0, 


where  Xj-ix,^)  positions  in  the  image,  or  output,  plane  of  the  imaging  system 
and  M  is  the  lateral  magnification  of  the  system.  The  intensity  in  the  image  plane 
is  the  coherent  superposition  of  the  wavefronts  B\  and  B2  such  that 

HB  ,+fl  xWBA’+Bfi'i+B  jB*,. 

As  always,  the  last  two  terms  represent  the  interference  of  the  two  wavefronts.  To 
separate  the  microscopic  (speckle)  interference  from  the  macroscopic  interference 
it  is  useful  to  examine  the  average  interference  field  over  a  large  number  of 
diffusers,  p.  The  emsemble  average  of  the  interference  terms  may  be  expressed  as 


<2?  xBi^cB  i  >. 


However,  since  <BxB\>  —  <B-^\  >  ,  the  overall  average  will  be  twice  the  aver¬ 
age  of  either  term.  Therefore  consider 


<B\B\  >— <( J'ai(xo)piOco)Mx)- Mx<i)dxo)x{f  a\(yo)p&£o)h' (x,— Mx0)dxo)>. 
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By  changing  the  dummy  variable  in  the  second  integral  to  x0',  the  order  of  multi¬ 
plication  and  integration  may  be  reversed  such  that 


<B\B\  >—<ffa  i(xo)a\(xQ)pi(  l))p2(xo  )h(jcJ—Mxc,'/t'(xJ—Mx0‘)dx^ixQ'>. 


Finally,  p\  and  p2  are  considered  n  be  ergodic  random  processes  which  is  to  say 
that  the  statistical  properties  of  th  roughness  over  any  sample  area  of  the  object 
surface  are  equal  to  those  obtained  for  any  other  sample  area.  Consequently  it  is 
permissible  to  reverse  the  order  v.  h  which  the  integration  and  averaging  is  per¬ 
formed  (Lathi,  1968)  so  that 

<B\B'i>— ffa  i(xo)fliO£o')<Pi(  )pl(x0')>h Cx,  —Mx0)h * (xJ—Mx0''ktxodxol  19) 


Through  Equation  19  it  is  possii  e  to  demonstrate  a  great  deal  about  the  effect 
of  surface  roughness  and  speckle  m  the  interference  of  two  diffusely  reflected 
wavefronts.  First,  to  convince  oneself  of  the  physical  significance  of  Equation  19, 
consider  the  case  where  B^Bi,  then 

<BXB\  >-<ByB\  >-<|5,|2>, 

which  is  the  average  intensity  of  B j  in  the  output  plane  of  the  imaging  system. 
Since  Bf-Bi  then  Pi”p2.  Thus  <Pi(xo)p*(xo)>  is  the  autocorrelation  of  the  sur¬ 
face  roughness  Rp{xo—xo')-  Recall  that  the  spatial  frequency  spectrum  of  the  sur¬ 
face  roughness  is  assumed  to  extend  well  beyond  the  bandwidth  of  the  imaging 
system.  In  other  words  the  fine  detail  of  the  surface  cannot  be  resolved  by  the  im- 


aging  system.  Therefore  the  autocorrelation  peak  of  the  surface  roughness  will  be 
much  narrower  than  the  resolution  of  the  imaging  system  and  can  be  approximated 
by  a  delta  function: 

<P\(xd)p\(xo')>=:s&(x(rxo)- 

Owing  to  the  shifting  property  of  the  delta  function,  Equation  19  is  considerably 
simplified, 

<|5|2>— ffa\(x0)a]  (xo)h(Xj-Mxo) 

x  ri '  (x;  —Mx0)&(xQ—x0')dxtfix0' 

-J"l  a  ]\1ix0)\h\Kx1-Mxo)dX(j. 

The  above  expression  is  now  recognized  as  the  convolution  of  the  input  intensity, 
|a|J,  with  the  incoherent  transfer  function. 

Having  established  the  physical  significance  of  Equation  19,  return  now  to  the 
interference  case  where  First  observe  that,  for  interference  to  exist  at  all, 

the  roughness  terms  for  the  two  wavefronts  must  be  correlated,  i.e. 
<Pi(xn)P7(xn')>g*0  •  In  practice,  then,  one  can  interfere  two  speckle  fields  only 
if  both  wavefronts  originate  by  reflection  from  the  same  object  surface  in  two 
different  states.  This  explains  the  value  of  holography  in  performing  inier- 
ferometnc  analysis  of  diffuse  object  surfaces.  Assume  that  following  the  recording 
and  subsequent  development  of  a  hologram,  the  test  object  is  translated  by  an 
amount  u  as  shown  in  the  real  time  holographic  interferometer  of  Figure  11.  The 
interference  terra  (Equation  19)  will  be  affected  in  two  ways.  First,  Bi  will  be 
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Figure  11:  Real  time  holographic  interferometer. 
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shifted  in  phase  by  an  amount  equal  to  the  projection  of  u  on  the  sensitivity  vec¬ 
tor,  K  (see  Figure  9).  Secondly,  there  will  be  a  lateral  shift  in  the  speckle  field  of 
B2  in  the  image  plane  by  an  amount  equal  to  the  product  of  the  lateral 
magnification,  M,  and  the  projection  of  u  onx<>  Thus 

where 

tl  —M (u  xqM^K  u. 

Substituting  into  Equation  18: 

Bf-ej4fa  \(xo)pi(xo)h  (x,  +U,  —Mx<ddx0. 

So  that  Equation  19  becomes 

<B\B\  J*a  i(x0)ai  (*o)  <  P^x.o)p'(xo)>h  (x,  —Mx0)h  *  (x,  -hr,  —Mx0)dxodxo 

^*SSa  \(x0)a\(xo)k(^-Mx0)h'(^+Uj~Mx0Mx0-xo)dxodx0 
Ha  il  Kxq  )eJ*f h  (x,  —Mx0  )h’  (x,+Uj—Mxo)dxo' 

~\o  \\^x0'kJ4Rfl(uJ ), 

where  R*(u,)  is  the  autocorellation  of  the  coherent  impulse  response,  h.  This  au¬ 
tocorrelation  is  also  the  impulse  response  function  or  point  spread  diameter  and 
determines  the  resolution  of  the  imaging  system  which  is  the  speckle  size.  For  a 


circular  aperture,  the  autocorrelation  function  is  an  Airy  function  whose  value  is  a 
maximum  for  and  fails  to  zero  when  u,—  l.22k/D.  Therefore  the  interfer¬ 
ence  fringes  will  be  visible  only  for  lateral  shifts  of  the  speckle  pattern,  u, ,  which 
are  less  than  the  resolution  limit  of  the  imaging  system.  The  fringe  contrast,  y, 
can  then  be  defined; 


(20) 

With  the  development  of  the  concept  of  fringe  contrast  it  is  now  possible  to 
address  the  effect  of  speckle  on  the  accuracy  of  phase  measurements  in  heterodyne 
holographic  interferometry.  Since  y  takes  into  account  effects  of  the  imaging  sys¬ 
tem,  it  is  no  longer  necessary  to  define  the  output  wavefronts  B j  and  B2  in  terms 
of  the  input  fields  and  the  system  transfer  function.  Instead  let 


Bi(x,ty~’b](,x)vi(x)eJU'1  (21) 

B2  (x  ,r )— b  2(x  )v  2(x)eJ<'*1‘  +*(x)\ 


where  £>i  and  b2  are  the  average  macroscopic  image  amplitudes  and  vt  and  V2  are 
complex  variables  related  to  the  object  surface  roughness  and  the  resulting  speckle 
fields.  In  order  for  B\  and  B2  to  interfere,  however  it  is  necessary  that  V]  and  v2 
be  somewhat  correlated,  thus  v2  can  be  considered  to  be  the  sum  of  two  uncorre- 
lated  fields  whose  magnitudes  are  dependent  on  the  contrast,  y  (Dandliker,  1980); 


v2(x  >-yv  \(x)eJ*+{  1— y2)*\v(x) 


<ViV^  >—/$»(«) 

<v,w*  >H) 
<|v,l4>-2. 


Substituting  into  Equation  21: 

B  |(x  ,t )— i b  fa  )v  (x  )ej" 

Bfe,t)-b1{x)[yv(x)eJ*M<\—yiy,w(jc)]eJmif. 

The  interference  term  of  the  intensity: 

B\B1a[yvv'eJ++{l-y1)*v'w\eJ<mi~m')‘,  (22) 

where  the  first  term  is  the  correlated  term  and  the  second  is  an  uncorrelated  term 
which  will  add  randomly  to  the  coherent  term  introducing  a  phase  error  For  a 
detector  area  large  enough  only  to  receive  light  from  two  speckles,  the  vector  di¬ 
agram  in  Figure  12  applies.  The  net  uncorrelated  component  is  the  vector  sum  of 
the  random  uncorrelated  components  of  each  speckle.  The  correlated  components 
add  coherently  resulting  in  a  net  component  equal  to  the  sum  of  the  individual 
components  from  each  speckle.  One  can  see  then  that  the  phase  error,  Ad>,  will  be 
reduced  as  the  number  of  speckles  falling  on  the  detector  increases.  From  the 
figure,  the  phase  error  can  be  expresses  as 


For  small  phase  errors  sinA4*~ A$,  so  the  mean  square  phase  error  is 


<&t2> 


<|/J3sinV>_i  <|/J2> 

<|/cl2>  2  <|/c|2> 


To  compute  A&  as  a  function  of  the  number  of  speckles,  N,  on  the  detector,  it  is 
necessary  to  find  expressions  for  <|jc|2>  and  <|/J2>.  From  Equation  22 

<l4l2>-  Z  <wi  > 

11— m— 1 

N 

-T2  Z  <V„vy„Vm> 

«— m— 1 

-y\N  <|v  |  >+JVOV-l )  <|v  |2>2 
-y2M!V+ 1); 

<l4l2>*“  Z  > 

it— m— 1 

V 

— {  1— y2)  Z  <vn”rVm*’m> 

n—m—\ 

-<  1— y2)JV<|v|2><|w|2> 

— ( l— y2)JV , 


and  finally: 


<A*2> 


i 


1  <l/ul2> 

2  <|/cl2> 

1  JV(  1-y2) 

2  y2JV(JV+l) 

(1-y2) 


I 


The  use  of  two  reference  beams  to  record  a  double  exposure  hologram  for  sub¬ 
sequent  heterodyne  analysis  introduces  two  additional  parameters  whose  effects 
may  degrade  the  accuracy  of  the  phase  measurements.  These  two  parameters  are 
misalignment  of  the  reconstructing  reference  beams  and  wavelength  differences 
between  the  recording  and  reconstruction  sources.  Both  are  pertinent  to  this  work 
since  in  the  final  experiments  holograms  were  recorded  and  analyzed  in  two 
separate  optical  systems.  In  addition  a  frequency-doubled  NdtYAG  laser 
(A— 532nm)  was  used  to  record  the  holograms  while  reconstruction  and  analysis 
was  performed  using  the  S14.5nm  line  of  an  argon  ion  laser. 

Although  Dandliker  discusses  effects  of  misalignment  and  wavelength  in  his 
1980  review  paper,  the  physical  significance  of  these  effects  is  easily  lost  in  the 
mathematics  associated  with  k-space  manipulations  and  scalar  diffraction  theory. 
For  this  reason,  the  following  development  is  offered  using  geometrical  optics 
theory  of  diffraction  gratings.  While  this  approach  does  not  appear  in  the  litera¬ 
ture,  it  is  never-the-less  appropriate  for  this  background  discussion. 

In  one  sense,  a  hologram  is  a  spatial-frequency  moduated  diffraction  grating 
whose  carrier  frequency  is  determined  by  the  angle  between  the  reference  beam 
and  the  average  object  beam  direction.  In  the  case  of  a  dual  reference  beam  holo- 


gam,  two  independent  diffraction  gratings  are  recorded  on  the  hologram  -  one  for 
each  exposure.  Given  that  the  hologram  may  be  treated  as  a  complex  diffraction 
gating,  it  becomes  clear  then  that  any  angular  misalignment  or  wavelength  change 
will  cause  a  change  in  the  angles  of  the  diffracted  beams  exiting  the  hologram.  An¬ 
gular  changes  result  in  lateral  shifts  of  the  images  formed  in  the  output  plane  of 
the  interferometer.  Such  lateral  shifts  cause  a  reduction  in  fringe  contrast  accord¬ 
ing  to  Equation  10  which  in  turn  degrades  the  phase  accuracy  as  shown  in  Equa¬ 
tion  23. 

To  quantitate  these  effects,  consider  the  hologram  geometry  shown  in  Figure 
13  where  R\  and  R2  are  plane  waves.  The  hologram  is  dawn  as  a  broken  line 
symbolic  of  its  diffraction  grating  nature.  An  imaging  lens  is  located  as  close  to 
the  hologram  plate  as  possible.  Assuming  perfect  alignment  of  the  reconstructiang 
reference  beams,  Rl  and  R2,  and  further  assuming  no  change  in  optical 
wavelength  between  hologram  recording  and  reconstruction,  the  reconstructed  ob¬ 
ject  beams,  0\  and  02,  propagate  at  an  angle  9j  to  the  hologram  (grating)  normal 
and  are  imaged  on  the  output  plane  where  they  interfere.  The  propagation  angles 
for  R  i,  O  i,  R2,  and  02  obey  the  grating  equation  for  first  order  diffraction; 

A-xr  (sin  0/j+sinftj )» 


where  the  0  angles  are  measured  with  respect  to  the  gating  normal  and  a  is  the 
gating  spacing.  In  the' case  diagramed  in  Figure  13; 


A— a  ,(s//i(0|— ' 8)+»n  9r) 


(24) 
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k—a  2(sjn  (  hhan  9r ) 

Consider  the  case  when  both  reference  beams  are  slightly  misaligned  such  that 
each  is  displaced  angularly  from  its  bisector  by  an  additional  amount  AS.  In  this 
case  for  R\  the  change  in  9r  may  be  computed  by  direct  application  of  Equation 
24: 


sinfy— sin(0/+-AS,  W  — — sn  (9— 6)]— { — — sin(0— (6+A8))] 
ai  o, 

Icos^  (20/4-Afy  )sin~i—^9i  >-sin(0— 5— AS )— sin(  9—6) 

-2cos^  20— 28— A8)smy(— AS). 


If  A8  is  small,  then 


—A 9r  c  ;s 9r  =— AScos  (  0—8  ) 


A  9 


AS  cos  (9—6) 

=“A8— - - - 

cos 9t 


Similarly  for  R2, 


A  9,- 


A  *  cos(9—6) 

AS - - - 

cos  9 1 


The  total  angular  separation  between  the  two  reconstructed  object  beams  is 


AS/  “AS/j — AO/j 

— A8 — -  [coa(0 — 8)+-cos(0+8)] 
cos  9, 


- ^-[2ms4(-28)cos^(2fl)] 

cos  9j  2  2 


-2A8 


cosfl 

COS^f 


x)s8 


— 2A8cos8  since  9r—6. 


The  corresponding  translation  in  the  image  plane  is 

U/—  Afyz 
—2zA8cos8. 


(25) 


For  a  practical  system  using  a  circular  imaging  aperture  so  that 


y 


2  jf™'0) 

7TUfD/k2  1  Az 


and  A— 514.5nm,  z—  15cm,  D— 7.5cm,  8— 5°,  the  maximum  angular  misalignment, 
AS,  which  will  still  result  in  a  fringe  contrast  y-0.7  may  be  computed  since  y— 0.7 
corresponds  to  ufZ)/Az-0.5.  By  substitution,  Uj  is  found  to  be  equal  to  the 
wavelength  (514.5nm).  Solving  Equation  25  for  AS  ,  the  maximum  angular 
misalignment  is  found  to  be  about  lCT^degreesl  Such  fine  alignments  are  achiev¬ 
able  with  precision  optical  mounts.  However,  this  result  helps  explain  why  fre¬ 
quent  realignment  is  required  with  changes  in  the  ambient  temperature  of  the  la¬ 
boratory. 


Sensitivity  to  wavelength  changes  can  likewise  be  analyzed  using  the  diffraction 
grating  analogy  and  Figure  13.  This  time  it  will  be  assumed  that  ASH)  and  that  a 
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different  wavelength  is  used  for  reconstructing  the  holographic  images  than  was 
used  to  record  the  hologram.  The  reconstruction  angle,  9r,  is  shifted  when  the 
hologram  is  reconstructed  using  X2  from  what  it  would  have  been  if  reconstructed 
using  Xj  by  an  amount  which  may  be  computed  using  Equation  24: 


sinfy— s//i(0/+A0/ )— [  — — am (g — 6)1 — C  ■  2  —an {9— 8)] 
o  i  a2 

— Afy  cos0/“ — (  X  i — X2) 


A0;  - ^-—-(X)— ^ X2). 

1  a  icos  9f 


M 

"M 


The  same  result  is  obtained  for  the  second  reference  beam  indicating  as  expected 
that  9r  will  be  shifted  in  the  same  direction  for  both  cases.  Therefore  the  total 
shift  in  0r  is  the  difference  between  the  two: 


A0,“A0;]-A0/2 

^1~^2  ^1~*2 
a  )Cos Oj  a  icos Of 
_  X|— X2  Q|— a2 
cos  0r  fl!a2 


From  Bragg’s  equation 


2s//i-(20-8) 


2s//iy(20-8) 

_ _  X,2»>i^(2W-«>-X12aj>i-J-(2e-«) 

a  i  —a  2  2  2 


«.M.a  i_  .  r.«»  «... x 


-OS,  f. — •  *--•»  ■-  mAj.1  «  A  m-  •"  M  ' 


.-.'.VaV-  . 


Pulsed  Holography 


Holographic  recording  using  pulsed  lasers  has  been  used  to  reduce  the  effects 
of  environmental  vibration  or  object  motion  on  holographic  quality  over  a  broad 
range  of  applications  from  holographic  portraiture  to  nondestructive  testing  (Col¬ 
lier,  1971b).  Nondestructive  testing  applications  include  contouring,  vibrational 
mode  analysis,  and  transient  stress  and  deformation  analysis  (Gagosz,  1974;  Chu, 
1972;  Albe,  1976).  In  the  area  of  transient  analysis,  pulsed  holography  has  been 
employed  to  "freeze"  the  motion  of  relatively  large  amplitude  acoustic  waves  in¬ 
cluding  plate  waves  and  Rayleigh  waves  (Apia ham ian,  1971;  Gagosz,  1974).  In 
both  of  these  studies,  however,  the  surface  displacements  were  sufficient  to  be 
visible  using  single  reference  beam  double  exposure  holographic  interferometry. 
They  are  significant  to  this  study  in  that  they  demonstrate  that  surface  displace¬ 
ments,  and  rate  of  displacement  in  particular,  associated  with  acoustic  waves  is 
insufficient  to  cause  noticable  reduction  of  holographic  efficiency.  Indeed  the 
reduction  of  holographic  efficiency  may  be  computed  directly  for  objects  in  uni¬ 
form  and/or  vibrational  motion  if  the  laser  pulse  duration,  pulse  shape,  and  the 
component  of  object  velocity  along  the  sensitivity  vector  is  known  (Gustafson, 
1980).  Using  Gustafson’s  formulas,  the  reduction  of  holographic  efficiency  result¬ 
ing  from  the  subfringe  displacements  encountered  in  this  work  would  be  less  than 


Although  the  reconstruction  efficiency  is  not  significantly  reduced  by  small 
acoustic  perturbations  of  the  object  surface,  the  accuracy  with  which  such  perturba¬ 
tions  are  recorded  remains  in  question.  Unfortunately  the  answer  to  this  question 
is  not  available  in  the  literature  at  present.  However  an  argument  can  be  present¬ 
ed  based  on  the  supposition  that  one  cannot  expect  recorded  phase  accuracy  any 
greater  than  that  associated  with  the  object  displacement  which  takes  place  during 
the  laser  pulse.  To  get  a  worst-case  estimate  of  phase  error  using  this  argument  it 
will  be  necessary  to  estimate  the  surface  velocity  normal  to  the  object  surface  as  an 
acoustic  wave  propagates  across  it.  For  the  sake  of  simplicity,  assume  that  the  nor¬ 
mal  (z)  displacement,  w(x,t)  of  the  object  surface  resulting  from  a  surface  wave 
propagating  in  the  x-direction  may  be  described  as  a  sinusoid  as  follows: 


w(x  ,t  )-—-bV(pin  (  an  —kgX  ), 


where  ka  is  the  acoustic  wave  propagation  constant  2tt/A*  and  o>  is  the  acoustic 
frequency.  The  displacement  velocity,  then  is: 

dw 

— -wqCOs  ( <ut—kax  ), 

at 

the  maximum  of  which  will  occur  when  Thus 

dw. 

-^iraax-^OW- 


An  acoustic  frequency  of  500kHz  corresponds  to  an  acoustic  wavelength  of  about 


6mm  for  aluminum  where  the  surface  wave  velocity  is  92%  of  the  bulk  shear  wave 
velocity.  In  the  experiments  documented  here,  a  frequency  doubled  NdtYAG 
laser  (X-*532nm)  with  a  half  amplitude  pulse  duration  of  9ns  is  used.  Therefore 
the  phase  uncertainty,  A<£,  associated  with  the  finite  recording  pulse  width  is 

2k{w0at)T  (27) 

"2s^r»(2<rS00x  lo3)(9* ur*)H'0 

— 6.68x  IOVq  radians 

*-3.83x  107w0  degrees. 


Thus  for  a  peak  wave  amplitude  wq— y  or  266nm,  the  phase  uncertainty  1°. 

Bear  in  mind  that  a  surface  displacement  as  large  as  266nm  will  produce  one  holo¬ 
graphic  fringe  in  the  interference  field.  Consequently  a  1°  phase  error  would  still 
permit  measurement  of  the  wave  peak  to  one  part  in  360  or  about  0.3%.  Accord¬ 
ing  to  Equation  27,  one  may  expect  better  base  accuracy  for  smaller  surface  dis¬ 


placements. 


Experimental  Design 


Three  primary  experiments  were  performed  to  demonstrate  that  heterodyne 
analysis  of  pulsed  holograms  can  be  used  to  measure  sub-fringe  displacements  as¬ 
sociated  with  traveling  surface  acoustic  waves.  In  the  first  experiment,  a  CW 
heterodyne  system  designed  very  much  after  Dandliker’s  was  built  Using  this  sys¬ 
tem,  the  processing  electronics  were  optimized  and  sub-fringe  interpolation  to 
1/900  of  a  fringe  (2.9a)  was  demonstrated.  The  last  two  experiments  addressed 
the  problems  of  synchronization  of  the  Laser  output  with  the  sound  wave  source 
for  multiple  exposure  holography.  As  part  of  the  third  and  final  experiment,  holo¬ 
grams  recorded  using  the  frequency  doubled  NdtYAG  laser  were  analyzd  using 
heterodyne  readout  and  phase  detection. 


CW  Heterodyne  Holographic  Interferometry 


The  holographic  recording  system  diagramed  in  Figure  14  was  used  to  record  a 
double  exposure,  dual  reference  beam  hologram  of  the  articulating  surfaces  of  the 
tibia!  component  of  an  artificial  knee  joint  (see  Figure  15).  The  two  articulating 
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Figure  14:  CV7  recording  set-up  for 
heterodyne  holography. 
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surfaces  are  concave  and  serve  as  sockets  into  which  the  condyles  of  the  upper, 
femoral,  component  are  located.  The  object  was  placed  in  a  sealed  chamber  such 
that  the  surfaces  of  interest  could  be  viewed  through  the  optically  flat  front  win¬ 
dow  (see  Figure  16).  A  contour  hologram  of  the  object  was  recorded  using  the 
multiple  refractive  index  technique  referenced  earlier.  %  changing  the  pressure 
within  the  test  chamber  between  exposures,  the  contour  interval  could  be  continu¬ 
ously  adjusted  so  that  less  then  one  contour  fringe  could  be  observed  over  the  im¬ 
age.  Thus  the  details  of  the  surface  contour  were  revealed  only  by  sub-fringe  inter¬ 
polation. 

When  recording  a  hologram,  the  beam  from  a  Spectra-Physics  165  argon  ion 
laser  was  split  and  directed  through  three  acousto-optic  modulators  (refer  to  Figure 
14).  Each  modulator  (Isomet  1205C)  was  driven  by  an  80 MHz  signal  from  a  sin¬ 
gle  oscillator/driver  during  exposure.  Beam  ratios  were  adjusted  so  that  either 
reference  beam  was  one  to  two  times  as  bright  as  the  object  beam  at  the  film 
plane.  A  different  reference  beam  was  blocked  for  each  exposure.  Agfa  10E56- 
N  AH  film  in  a  4"x5"  glass  plate  format  was  used.  Total  exposure  was  approximate¬ 
ly  10  ergs/cm  2.  The  film  was  developed  for  5  minutes  at  room  temperature  using 
Kodak  EM 9  developer.  After  a  30  second  rinse,  the  plate  was  fixed  with  Kodak 
rapid  fixer  for  5  minutes. 

Once  dty,  the  film  plate  was  replaced  in  its  holder  for  reconstruction  and 
analysis.  The  beamsplitter  used  to  provide  illumination  of  the  object  was  removed 
as  was  the  etaion  within  the  laser  cavity.  Both  these  steps  increased  the  amount  of 
light  available  in  the  reference  beams.  Precise  replacement  of  the  film  plate  was  ac- 
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Figure  16:  Test  charier  for  contouring 
knee  implant. 
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corapiished  by  observing  the  interference  pattern  in  the  film  plane  and  adjusting 
the  film  position  and  final  reference  mirrors  to  minimize  the  fringe  order.  Once  the 
film  plate  was  repositioned,  the  drive  signal  from  one  of  the  acousto-optic  modula¬ 
tors  was  changed  from  the  common  80 MHz  oscillator  to  a  second  oscillator  this 
time  running  at  80.1MHz.  Merit  model  R1S00  silicon  PTN  diodes  with  integral 
transconductance  amplifiers  were  used  for  signal  detection.  The  reference  detector 
was  positioned  for  best  S/N  ratio.  A  second  detector  was  mounted  on  an  x-y 
translation  stage  driven  by  Newport  Research  model  850-1  linear  actuators.  The 
reference  and  probe  signals  were  amplified  by  Tektronix  AM501  differential 
amplifiers,  passively  band-pass  filtered  at  100 kHz  and  fed  to  the  reference  and  sig¬ 
nal  inputs  of  a  PAR  model  5301  lock-in  amplifier.  The  lock-in  amp  was  configured 
to  measure  phase  difference  between  the  reference  and  probe  signals.  The  analog 
phase  output  signal  was  connected  to  the  Y-axis  input  of  a  Hewlett-Packard  7045  A 
x-y  recorder.  The  x-axis  was  driven  with  a  ramp  signal  from  an  HP  3304  function 
generator.  (In  subsequent  experiments  probe  detector  positioning,  phase  voltage 
recording,  processing,  and  display  were  ail  performed  under  computer  control.  The 
implementation  of  this  computerized  data  acquisition  and  analysis  system  will  be 
discussed  later.)  Using  the  actuator  joystick  controller,  the  probe  detector  was  po¬ 
sitioned  to  one  side  of  the  output  image  and  then  scanned  across  the  image  at  a 
constant  velocity.  The  beginning  of  the  probe  scan  was  synchronized  with  the 
recorders  x-axis  drive  signal.  Ry  incrementing  both  the  recorder  pen  and  the 
detector  probe  along  the  y-axis,  a  y-modulated  raster  scan  of  the  knee  prosthesis 


was  generated. 
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Double  Pulsed,  Single  Reference  Beam  Holography 


A  Quantei  model  580-C  NcLYAi}  Laser  with  an  etaion  and  frequency  doubling 
crystal  was  used  to  record  all  of  the  pulsed  holograms.  The  laser  has  a  horizontally 
polarized  output  with  a  pulse  width  ;f  9ns  and  an  output  energy  of  1.2  joules/pulse 
at  1064nm  and  about  half  that  energy  at  532nm.  The  laser  energy  in  the  green 
(53* Tm)  can  be  reduced  either  by  tentional  misalignment  of  the  doubling  crystal 
or  by  disabling  the  amplifier  flash  mps.  In  the  later  case,  the  maximum  energy 
per  pulse  at  S32nm  from  the  osc  uor  alone  is  350mJ.  Tt  was  this  low  evergy 
configuration  which  was  used  to  n  ird  holograms  of  specularly  reflecting  objects. 
Recording  of  diffuse  object  surface  oquired  higher  illumination  intensities. 

The  recording  geometry  diagrar  <1  in  Figure  17  was  used  to  record  the  travel¬ 
ing  surface  wavefront  propagating  ac  -ss  an  aluminum  plate.  The  plate  dimensions 
were  25cm  x  38cm  x  3.8cm  thick.  Trie  front  surface  of  the  plate  was  grit  blasted 
to  provide  a  diffusely  reflecting  surface  The  glass  plate  beamsplitter  was  posi¬ 
tioned  near  the  brewster  angle  to  permit  most  of  the  beam  energy  to  pass  through 
it  to  illuminate  the  object. 

The  purpose  of  this  experiment  was  in  part  to  verify  the  timing  and  allignment 
of  the  pulsed  holographic  recording  system.  Since  no  fringe  interpolation  was  to 
be  done  at  this  stage  it  was  necessary  that  the  recorded  acoustic  wavefronts  be 


sufficiently  large  to  produce  interference  fringes  which  may  be  visually  observed 
when  the  hologram  is  reconstructed.  In  order  for  this  to  occur  the  the  surface  wave 
displacement  amplitude  must  be  greater  than  A/2  or  about  2500 A.  To  generate 
such  large  surface  displacements  a  small  explosive  charge  was  used  as  the  sound 
source.  The  explosive  used  was  silver  acetylide  which  is  a  white  precipitate  formed 
when  acetylene  gas  is  bubbled  through  an  aqueous  solution  of  silver  nitrate.  Once 
rinsed  and  dried,  the  silver  acetylide  powder  was  pressed  into  small  pellets  about 
3mm  in  diameter  and  0.5mm  thick.  Each  pellet  was  then  sandwiched  between 
layers  of  cellophane  tape  and  taped  to  the  surface  of  the  plate  at  the  leftmost  edge. 
The  infrared  output  of  the  pulsed  laser  was  focused  on  the  pellet  to  detonate  the 
explosion. 

Assuming  that  the  surface  wave  velocity  is  about  92%  of  the  shear  wave  veloci¬ 
ty  in  aluminum  (3100  m/s),  then  the  surface  wavefront  takes  about  45/xs  to  pro¬ 
pagate  half  way  across  the  plate  surface.  Since  the  same  laser  was  used  to  detonate 
the  explosion  and  record  the  hologram,  it  was  necessary  to  modify  the  pock  els  cell 
triggering  circuit  to  double  pulse  the  laser  during  a  single  flash  lamp  pulse  with  a 
45/xs  pulse  separation.  The  pockels  cell  trigger  circuit  (Figure  18)  uses  eleven 
transistors  in  series  which  operate  in  avalanche  mode  to  deliver  a  high  voltage 
from  capacitor  C]  to  the  pockels  cell.  Once  avalanche  is  initiated  the  voltage 
across  the  series  transistors  must  fall  well  below  the  avalanche  threshold  voltage 
before  they  will  again  turn  off.  It  is  the  time  required  for  this  voltage  decay  which 
limits  the  rate  at  which  the  pockels  cell  may  be  triggered.  To  decrease  the  time  re¬ 
quired  to  reset  the  series  transistors,  an  additional  transistor  was  inserted  across 
the  pockels  cell  output  terminals.  This  transistor  turns  on  shortly  after  the  high 
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voltage  is  applied  to  the  cell  Once  on,  it  provides  a  low  resistance  path  which 
more  rapidly  bleeds  charge  from  the  capacitor  thus  causing  a  corresponding  in¬ 
crease  in  the  voltage  decay  rate.  Once  implemented,  this  modification  permits 
double  pulsing  of  the  laser  with  pulse  separations  externally  controlled  over  a 
range  from  20/us  to  100ms. 

Optimum  double  pulsing  of  the  laser  required  proper  liming  of  the  Q-switch 
pulses  relative  to  the  beginning  of  the  flashlamp  pulse.  For  this  reason,  the  trigger 
system  shown  in  Figure  19  was  used.  An  adjustable  digital  delay  generator  (BNC 
model  7010)  was  used  to  trigger  an  HP  model  8005 R  pulse  generator  with  double 
pulse  output.  The  delay  unit  was  triggered  by  the  flashlamp  drive  signal  generated 
internally  by  the  laser  control  electronics.  A  simple  buffer  circuit  capable  of  provid¬ 
ing  the  required  IS  volt  Q-switch  trigger  pulses  into  the  SO  ohm  trigger  input  cir¬ 
cuit  was  designed  and  built  since  the  TTL  output  of  the  HP  pulse  generator  was 
insufficient  to  trigger  the  Q-switch. 

The  sequence  for  recording  the  pulsed  acoustic  surface  wavefront  was  relatively 
straight  forward.  The  pulse  separation  was  set  for  45ms  and  the  digital  delay  adjust¬ 
ed  to  provide  equal  intensities  in  both  pulses.  The  explosive  charge  was  fixed  to 
the  plate  surface  and  a  film  plate  placed  in  the  holder.  When  the  laser  was  fired,  a 
holographic  exposure  of  the  plate  surface  was  recorded  using  the  green  line  of  the 
laser  output.  At  the  same  time,  the  light  from  the  infrared  port  of  the  laser  de¬ 
tonated  the  explosive  charge.  Forty-five  microseconds  later,  the  laser  emitted  a 
second  pulse  of  light.  This  time  a  second  exposure  of  the  hologram  was  made  us¬ 
ing  the  green  wavelength  so  that  any  surface  deformations  caused  by  the  explosion 
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Figure  19:  Laser  trigger  system 
and  timing. 


were  recorded.  The  infrared  beam  again  struck  the  plate.  However,  it  was  of  no 
consequence  since  this  time  there  was  no  charge  to  be  detonated.  Again  Agfa 
10E56  film  was  used  for  recording.  Kodak  D-19  developer  for  5  minutes  was  fol¬ 
lowed  by  a  water  rinse.  Instead  of  fixing  these  plates,  however,  the  plates  were 
beached  for  about  2  minutes  in  a  solution  of  potassium  dichromate  and  sulfuric 
add  to  improve  diffraction  efficiency. 


Surface  Acoustic  Wavefront  Recording  for  Heterodyne  Analysis 


Through  this  experiment,  the  attempt  to  measure  and  display  subfringe  detail 
of  displacements  associated  with  surface  wave  propagation  was  made.  The  optical 
arrangement  diagrammed  in  Figure  20  was  used  to  record  dual  reference  holo¬ 
grams  of  the  surface  wavefront  generated  by  detonation  of  a  silver  acetylide  pellet 
as  described  in  the  previous  experiment.  The  same  aluminum  spedmen  used  in 
the  previous  experiment  was  again  used  as  the  test  object.  This  time  however,  the 
object  surface  was  polished  to  reduce  scattering  of  the  collimated  illuminating 
beam. 

To  fadlitate  heterodyne  analysis,  a  dual  reference  beam  double  exposure  holo¬ 
gram  could  have  been  recorded  as  described  previously  in  the  Background  section. 
Using  a  double-pulsed  laser  to  both  detonate  the  expolsive  charge  and  expose  the 


-82- 


holograms,  however,  would  have  required  some  type  of  electro-optic  switch  to 
change  the  reference  beam  path  during  the  45*ts  delay  between  the  laser  pulses. 
Unfortunately,  such  a  switching  device  was  not  available  for  this  experiment.  In¬ 
stead,  a  dual  reference  beam,  triple  exposure  hologram  was  recorded.  The  first  ex¬ 
posure  was  made  using  reference  path  1  as  shown  in  Figure  20.  During  this  first 
exposure,  the  IR  port  of  the  laser  was  dosed  to  prevent  detonation  of  the  explo¬ 
sive  charge.  After  the  diverting  prism  was  removed,  the  IR  port  was  again  opened 
and  the  laser  was  double-pulsed  to  record  exposures  2  and  3  using  reference  path 
2. 


While  the  triple  exposure  technique  permits  heterodyne  analysis  of  the  holo¬ 
graphic  interferogram,  it  has  a  pronounced  effect  on  the  calibration  and  sensitivity 
of  the  analysis.  This  effect  can  he  understood  by  examining  the  mathematical 
description  of  the  output  of  the  tr  ale  exposure  dual  reference  hologram.  The  in¬ 
tensity  distribution  in  the  output  j.  one  of  the  interferometer  is  nearly  identical  to 
that  described  by  Equation  17  exce  jt  that  an  additional  image  field,  Oj,  must  also 
be  added  in.  Thus  the  intensity  in  the  output  plane  becomes 
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It  is  clear  from  the  final  form  of  this  expression,  that  the  measured  phase 
differences  using,  the  triple  exposure  technique  will  be  half  that  which  would  be 
measured  using  the  double  exposure  technique.  Further,  the  modulation  index 
will  be  zero  when  Atf  has  a  value  equal  to  odd  multiples  of  180°  and  will  be  nega¬ 
tive  for  1800<A^<540°  the  effect  of  this  negative  modulation  index  is  to  generate 
a  180°  discontinuity  in  the  output  of  the  phase  detector  when  the  optical  phase 
difference  between  corresponding  points  on  o  j  and  exceeds  180°.  Fortunately, 
for  subfringe  displacements,  optical  phase  differences  will  always  lie  between 
-t-180"  and  -180°  thus  insuring  that  the  modulation  index  will  always  be  positive. 
As  the  modulation  index  falls  below  1,  the  signal  to  noise  ratio  degrades  as 
descnbed  earlier.  Again  for  small  displacements,  this  may  not  be  a  real  problem 
since  M  will  always  be  between  0.8  and  1.0  for  optical  phase  differences  between 
0°  and  ±120°. 

Talcing  the  above  into  consideration  it  was  felt  that  the  triple  exposure  tech¬ 
nique  could  be  used  to  demonstrate  the  application  of  heterodyne  analysis  to 


pulsed  holographic  interferometry.  The  triple  exposure  hologram  was  recorded  on 
Agfa  10E56NAH  plates  and  processed  as  before  The  processed  plates  were  then 
placed  in  the  heterodyne  readout  system  diagrammed  in  Figure  21.  The  plates 
were  reversed  in  the  film  holder  so  that  the  conjugate,  real  image  was  projected  on 
the  output  plane  without  requiring  imaging  lenses.  Alignment  was  performed  with 
both  A-O  modulators  driven  at  80.0MHz  so  that  a  stationary  fringe  pattern  could 
be  observed.  Once  the  system  was  aligned,  the  modulators  were  driven  at 
80.0 MHz  and  80.1MHz,  respectively. 

Figure  22  shows  the  data  collection  system  in  block  diagram  form.  Again,  two 
Merit  R1500  optical  detectors  were  used  to  pick  up  the  100kHz  beat  signal.  This 
time,  however,  no  beamsplitter  was  used  to  provide  a  reference  signal.  Instead, 
the  reference  detector  was  placed  in  the  output  plane  at  one  corner  of  the  project¬ 
ed  image.  The  scanning  detector  was  mounted  on  an  x-y  translation  stage  driven 
by  two  NRC  85C-C  motorized  micrometers.  The  micrometers  were  driven  through 
an  interface  controlled  by  a  Rockwell  AIM  65  microcomputer.  All  micrometer 
drive  routines  were  written  using  FORTH  programming  language  and  are  listed  in 
the  Appendix.  After  selecting  the  appropriate  vertical  position  with  the  y-axis  mi¬ 
crometer,  the  horizontal  scan  routine  was  run.  This  routine  scanned  the  photo¬ 
detector  over  a  distance  of  about  1cm  in  125  Mm  increments.  After  each  increment, 
a  trigger  pulse  was  sent  by  the  AIM  65  to  the  A/D  converter  module  of  an  IBM 
model  7406-1  Device  Coupler.  The  trigger  pulse  signaled  the  module  to  convert 
the  analog  phase  difference  signal  from  the  lock-in  amp  and  to  transmit  that  digi¬ 
tized  value  through  an  RS232  interface  to  an  IBM  5100  computer.  The  5100  com¬ 
puter  was  used  to  decode  the  digitized  data  and  store  the  data  in  vector  form  such 
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that  each  90  element  vector  corresponded  to  one  scan  across  the  output  field.  The 
raw  data  received  by  the  computer  is  modulo  360°  in  that  phase  values  measured 
by  the  lock-in  amp  are  restricted  to  those  between  -180°  and  +180°.  Consequent¬ 
ly,  a  computer  routine  was  used  to  ‘unwrap"  the  raw  data  to  remove  these  large 
discontinuities.  As  predicted  wuh  the  relatively  large  surface  displacements  associ¬ 
ated  with  an  explosive  sound  source,  180°  discontinuities  in  the  data  are  also  ob¬ 
served.  These  too  are  removed  by  computer  processing.  Finally  the  processed 
data  was  plotted  on  a  Tektronix  4013  graphics  display  terminal  and  subsequently 
ploited  on  a  Hewlett-Packard  7045  A  x-y  recorder  using  the  Analog  Output  Module 
of  the  Device  Coupler.  All  collection  and  processing  software  in  the  IBM  5100 
computer  was  written  in  APL  and  is  listed  in  the  Appendix. 
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Results  and  Discussion 


Sub-fringe  Interpolation:  Contour  Holography 


The  image  shown  in  Figure  23  is  the  output  of  the  dual  reference  contour 
hologram  of  the  knee  implant  recorded  as  described  in  the  previous  chapter.  In 
order  to  photographically  record  the  contour  fringes  in  the  output  field,  both  A-O 
cells  in  the  reference  legs  were  driven  at  the  same  80.0MHz  drive  frequency.  Dur¬ 
ing  heterodyne  analysis,  the  two  modulating  frequencies  differ  by  100kHz  thus 
causing  the  fringe  pattern  to  fluctuate  at  such  a  high  rate  that  the  fringes  are  no 
longer  visible.  Optical  detectors  in  the  output  plane,  however,  detect  the 
sinusoidally  varying  intensity  and  transduce  it  into  a  voltage  signal  analog.  Initial¬ 
ly,  the  position  of  the  scanning  detector  relative  to  the  reference  detector  in  the 
output  plane  was  fixed  and  the  phase  difference  as  a  function  of  time  was  plotted 
in  order  to  measure  system  noise.  Using  this  technique,  the  short  term  phase 
noise  was  characterized  as  having  peak-io-peak  fluctuations  of  0.288°.  Long  term 
drift  of  about  0.04°  per  minute  was  also  observed.  In  subsequent  experiments,  the 
drift  rate  was  observed  to  deviate  from  this  value  indicating  that  drift  may  be  in 
part  the  result  of  dimensional  changes  in  the  apparatus  in  response  to  changes  in 


Figure  23:  Output  image  fo  Knee  implant 
when  =  w 2  =  80MHz. 


environmental  temperature. 


The  entrance  pupil  of  the  imaging  optics  was  the  film  plate  itself.  Therefore 

* 

I  the  speckle  size  was  determined  by  the  4  inch  ( 10cm)  dimension  of  the  film  plaie. 

Given  that  the  image  plane  was  some  50cm  behind  the  film  plate,  the  speckle  di¬ 
ameter  was  about  2.6  microns.  The  number  of  speckles,  N,  within  the  1mm  diam- 

I  eter  detector  was  then  about  1.5x  105.  From  Equation  23,  the  expected  phase  error 

with  iV— l.SxlO-5  and  an  estimated  modulation  index  of  0.8  was  about  0.11°. 
Therefore  the  overall  phase  noise  may  be  as  large  as  the  sum  of  the  elctronic  and 
speckle  noise  or  0.4°.  Tt  should  be  noted  that  this  0.4°  noise  level  corresponds  to 
an  optical  path  length  sensitivity  (for  SNR  —  1)  of  0.57nm  (5. 7 A)  or  an  object 
displacement  as  small  as  0.29nm  (2.9 A).  In  the  specific  case  of  contour  hologra¬ 
phy,  the  sensitivity  is  1/900  of  the  contour  fringe  interval. 

A  single  horizontal  line  scan  was  made  scanning  from  left  to  right  across  the 
image  plane  and  passing  through  what  would  be  the  center  of  the  stationary  fringe 
pattern  shown  in  Figure  23.  Note  that  the  plot  of  relative  phase  difference  as  a 
frunction  of  position  is  modulo  360°  (Figure  24)  and  that  if  "unwrapped"  it  would 
provide  a  profile  of  the  curved  surface.  However,  sufficient  sensitivity  existed  us¬ 
ing  this  analysis  system  to  detect  the  surface  features  of  interest  e'^en  when  the 
holographic  contour  interval  was  larger  than  the  depth  of  the  surface  relief  over 
the  entire  object.  Therefore,  a  second  contour  hologram  was  prepared  using  the 
double  refractive  index  technique  where,  as  described  previously,  only  a  slight 
change  in  air  pressure  inside  the  test  chamber  was  made  between  exposures.  The 
resulting  contour  image  showed  less  than  one  complete  fringe  over  the  object. 
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Thus  when  over  the  image,  the  phase  difference  between  the  scanning 

anl  refererae  signals  never  exceeded  360*.  Figure  23a  shows  a  series  of  horizontal 
scans  over  a  section  of  the  image.  The  spacing  between  the  scans  corresponds  to  a 
vertical  displacement  of  about  0.5mm  along  the  image.  The  scans  follow  the  nor¬ 
mal  curvature  of  the  artificial  knee  socket.  Just  short  of  halfway  across,  all  but  the 
lowest  scans  show  a  deviation  from  the  smooth  contour.  This  deviation  was 
caused  by  the  gouge  placed  in  the  socket  which  was  not  apparent  in  the  conven¬ 
tional  contour  hologram  (Figure  23).  A  section  of  the  gouge  is  shown  in  greater 
detail  in  Figure  25b  by  i rereading  the  gain  of  the  chart  recorder.  The  displacement 
between  scans  is  0.1mm  in  this  figure. 


Double  Pulsed,  Single  Reference  Beam  Holography 


Photos  of  the?  reconstructed  images  of  the  aluminum  plate  recorded  as 
described  in  the  corresponding  section  of  the  previous  chapter  are  shown  in  Fig¬ 
ures  26  a,  b,  and  c.  Figure  26a  shows  the  acoustic  wavefront  as  it  appeared  35  mi¬ 
croseconds  after  the  explosive  charge  was  detonated.  Figures  26b  and  26c  similarly 
show  wavefronts  from  subsequent  experiments  where  the  laser  pulse  separation 
was  40  and  45  microseconds  respectively. 


From  these  photographs,  the  wave  velocity  and  therfore  the  wave  type  can  be 
determined.  Table  T  summarizes  the  velocity  data  obtained  from  Figures  26  thus 
confirming  that  the  observed  primary  wavefronts  travel  at  a  velocity  corresponding 
to  the  Rayleigh  wave  velocity.  Following  the  primary  wavefront,  one  observes  an 
irregular  pattern  of  fringes.  This  pattern  varies  from  experiment  to  experiment  as 
the  laser  pulse  interval  is  varied.  Presumably  this  pattern  is  the  result  of  surface 
displacements  caused  by  the  superposition  of  waves  which  are  multiply  reflected 
from  the  front  and  bock  surfaces  of  the  test  plate. 


Heterodyne  Analysis  of  Surface  Acoustic  Wavefronts 


From  the  previous  experiment  it  was  confirmed  that  the  primary  surface  wave- 
front  passes  across  the  center  of  the  test  plate  some  45  microseconds  after  the  ex¬ 
plosive  charge  is  detonated.  This  single  data  point  was  important  since  those  holo¬ 
grams  recorded  for  subsequent  heterodyne  analysis  were  made  using  the  polished 
surface  of  the  specimen  as  discussed  earlier  and  only  an  8cm  diameter  region  was 
illuminated  by  the  object  beam,  ft  was  important,  therefore,  that  the  wavefront  be 
passing  through  the  8cm  circle  when  the  final  exposure  of  the  hologram  was  made. 


Two  steps  were  involved  in  the  alignment  of  the  developed  holograms  in  the 
heterodyne  analysis  system.  Driving  both  A-O  cells  at  the  same  frequency,  a  video 


35uS 

45iun 

2.66  x  10  cm/sec 

40uS 

54mn 

2.74  x  10^an/sec 

45uS 

59nm 

2.61  x  10^cm/sec 

Expected  Surface  Wave  Velocity  =  92%  x  3.1  x  105cm/sec  -  2. 85cm/ sec 
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camera  was  focussed  on  the  image  plane.  A  cross  which  marked  the  center  of  the 
test  surface  was  visible  in  the  images  reconstructed  by  both  reference  beams. 
Coarse  alignment  of  the  two  images  was  obtained  by  changing  the  two  reference 
beam  angles  to  superimpose  the  two  crosses  observed  on  the  video  monitor.  Final 
alignment  was  made  by  focussing  the  camera  on  the  plane  of  the  hologram  plate 
and  adjusting  the  reference  beam  angles  to  minimize  the  order  of  the  observed 
fringe  pattern.  As  predicted,  the  10s  separation  of  the  reconstructing  beams  and 
the  use  of  different  recording  and  playback  wavelengths  made  it  impossible  to  el¬ 
iminate  all  fringes  observed  in  the  hologram  plane.  Under  the  best  alignment  con¬ 
ditions,  less  than  one  fringe  was  observed  over  an  area  about  3cm  in  diameter  in 
the  center  of  the  illuminated,  region.  Beyond  this  central  region,  multiple  fringes 
encircling  the  center  were  observed.  Focussing  back  on  the  image  plane,  the  effect 
of  wavefront  misalignment  manifested  itself  in  reduced  fringe  contrast  and  fringe 
distortion  beyond  the  3cm  central  region.  Photographs  of  the  output  plane  with 
proper  and  improper  beam  alignment  are  shown  in  Figure  27.  Note  the  distortion 
of  the  fringes  in  Figure  27b  resulting  from  slight  misalignment. 

Figure  28a  shows  relative  phase  difference  as  a  function  of  position  across  a 
lcm  length  of  the  image.  This  raw  data  is  once  again  modulo  360*  and  exhibits 
the  180°  discontinuities  predicted  for  triple  exposure  heterodyne  holography.  Ac¬ 
counting  for  the  discontinuities  and  "unwrapping”  the  data  produces  the  plot  of  dis¬ 
placement  versus  position  of  Figure  28b.  Two  relatively  narrow  peaks  are  ob¬ 
served  in  the  leading  portion  of  this  wavefront  plot.  Although  it  was  beyond  the 
scope  of  this  research  to  explain  these  peaks,  it  was  originally  thought  that  they 
may  be  the  result  of  the  superposition  of  the  surface  wave  emanating  from  the 
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source  location  with  the  wavefront  which  is  reflected  from  the  edge  of  the  test 
plate  immediately  to  the  left  of  the  charge.  To  test  this  hypothesis,  a  second  ex¬ 
periment  was  performed  identical  with  the  first  except  that  the  charge  location  was 
displaced  towards  the  center  of  the  plate  a  full  centimeter  from  the  plate  edge.  As¬ 
suming  that  the  twin  peaks  were  the  result  of  interference  between  the  direct  and 
reflected  surface  waves,  the  wave  shape  and  especially  the  separation  between  the 
peaks  should  have  changed.  Instead,  as  can  be  seen  from  Figure  29,  both  the 
peaks  and  their  relative  spacing  remained  the  same  even  after  the  charge  location 
was  changed.  Therefore,  edge  reflection  does  not  account  for  the  dual  peaks  ob¬ 
served  in  the  surface  wave  shapes.  By  changing  vertical  positions  between  scans,  a 
full  field  taster  over  an  area  of  the  test  surface  can  be  obtained  and  plotted  in  "3- 
ET  as  shown  in  Figure  30. 

Several  sources  of  inaccuracy  limited  the  resolution  of  this  heterodyne  analysis 
system.  First,  the  measured  efficiency  of  the  triple  exposure  holograms  was  less 
than  1%  providing  an  output  field  intensity  of  about  1  microwatt/cm  .  From 
Equation  12,  the  expected  phase  variation  was  0.06s  for  a  fringe  contrast  of  0.7,  a 
lms  measurement  interval,  and  a  10kHz  bandwidth.  In  fact,  the  measured  signal 
variation  for  a  stationary  scanning  detector  was  0.6s  ,  which  was  much  worse  than 
predicted.  The  low  signal  intensity  in  the  output  plane  and  the  relatively  poor 
noise  figure  of  the  amplifiers  used  would  tend  to  indicate  that  the  detection  system 
was  not  shot  noise  limited  (a  key  assumption  in  deriving  Equation  12).  In  addition 
some  of  the  noise  may  be  attributed  to  variation  in  pathlength  owing  to  environ¬ 


mental  effects. 


A  second  noise  source  resulted  from  the  effort  to  conserve  light  in  the  output 
plane  by  eliminating  the  need  for  focussing  lenses.  As  a  result,  the  distance,  z, 
between  the  hologram  and  the  output  plane  was  a  fairly  long  61cm.  The  effective 
aperture  of  the  imaging  system  was  the  8cm  diameter  exposed  region  of  the  holo¬ 
gram.  The  speckle  size  was  determined  to  be  about  4.8  microns  in  diameter  so 
that  4.4x10*  speckles  were  encompassed  by  the  1mm  diameter  detector.  Using 
Equation  23  the  phase  error  resulting  from  speckle  noise  was  on  the  order  of 
0.28°. 

Finally,  the  rather  large  diameter  of  the  detector  itself  caused  some  spatial 
averaging  or  smoothing  of  the  surface  features.  As  a  result,  one  can  not  be  certain 
of  the  absolute  height  of  the  leading  edge  peaks  observed  in  Figures  31,  32,  and 
33.  It  is  quite  possible  that  these  pinks  are  considerably  larger  than  indicated. 
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Conclusions 


Subfringe  detail  of  traveling  surface  acoustic  waves  may  be  measured  and 
displayed  in  full  field  using  heterodyne  holographic  analysis  of  multiple  exposure 
holograms  recorded  with  a  pulsed  laser.  Systems  for  recording  and  displaying  sur¬ 
face  wavefronts  have  been  demonstrated  to  provide  displacement  resolution  to 
0.43nm  (4.3A)  over  a  3cm  diameter  field.  The  limited  field  size  is  the  result  of 
using  different  recording  and  playback  wavelengths.  The  field  size  may  be  im¬ 
proved  to  the  size  of  the  entire  field  recorded  by  the  hologram  by  changing  the 
recording  and  playback  geometries  to  make  the  two  reference  beams  more  nearly 
colinear.  Indeed  when  the  same  wavelength  was  used  for  both  recording  and 
analysis  of  the  contour  hologram,  analysis  of  the  entire  field  was  possible. 

The  single  most  important  step  m  improving  the  sensitivity  of  the  heterodyne 
analysis  system  would  be  to  increase  the  intensity  of  light  in  the  output  field. 
Several  options  are  available  to  increase  output  intensity.  First,  holographic 
efficiency  may  be  greatly  improved  by  using  a  different  recording  medium  either 
during  the  recording  of  the  hologram  itself  or  by  copying  the  primary  hologram 
onto  a  second  medium  using  a  longer  CW  copying  process.  Dichroma  ted  gelatins 
and  photoresist  holograms  have  been  reported  to  have  efficiencies  over  90%.  Using 
an  electro-optic  switch,  as  described  earlier,  to  allow  double  rather  than  triple  ex¬ 
posure  holograms  to  be  recorded  for  heterodyne  analysis  would  not  only  improve 
holographic  efficiency  but  would  provide  a  constant  fringe  visibility  in  the  output 


plane  independent  of  optical  phase  differences.  Still  more  improvement  could  be 
expected  if  a  higher  frequency  detector  capable  of  sensitive  recording  at  up  to 
80MHz  were  used.  This  would  permit  the  use  of  only  a  single  A-O  cell  for  fre¬ 
quency  shifting.  Using  the  undiffracted  beam  for  one  reference,  and  the  diffracted 
beam  for  the  other  reference,  virtually  all  of  the  light  from  the  laser  would  be  pro¬ 
jected  on  the  hologram.  Presently,  the  undiffracted  beams  from  the  two  A-O 
modulators  are  blocked  off  and  that  energy  is  lost. 

Increasing  the  light  available  in  the  output  plane  would  also  make  it  possible  to 
further  reduce  spatial  noise  in  the  system.  With  more  light  available,  a  lens  system 
may  be  used  to  reduce  the  working  distance  behind  the  hologram  and  thereby 
further  decrease  speckle  size  and  thus  improve  phase  accuracy.  In  addition,  the 
detector  area  may  be  reduced  to  further  improve  lateral  resolution  as  well. 

The  results  of  this  work  appear  to  support  hypotheses  proposed  by  the  author 
and  others  working  in  the  field  of  heterodyne  holographic  interferometry.  Rased 
on  the  experimental  results  obtained  and  the  mathematical  relationships  which 
have  been  developed  to  describe  system  performance,  the  author  feels  confident  in 
predicting  that  a  practical  system  for  studying  sub- Angstrom  transient  phenomena 
may  be  realized  using  heterodyne  holographic  interferometry. 
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